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Introduction: 

The  Laser-Materials  Interaction  Laboratory  at  the  Center  for  High 
Technology  Materials  of  the  University  of  New  Mexico  is  devoted  to  the  study 
of  a  broad  range  of  laser  spectroscopic  probes  of  semiconductor  and  nonlinear 
materials,  fabrication  processes  and  optoelectronic  devices.  Much  of  this 
work  is  being  carried  out  in  coiyunction  with  the  Optoelectronics  Research 
Center  Program  at  CHTM,  which  is  also  partially  funded  by  the  Air  Force 
Office  of  Scientific  Research.  Significant  progress  has  been  made  during  this 
reporting  period  in  a  number  of  areas  including:  ultrafast  operation  of 
optically-pumped  resonant  periodic-gain  GaAs  surface-emitting  lasers; 
grating  coupling  to  surface-plasma  waves;  and  high-speed  Si  Schottky 
photodiodes.  Only  brief  summaries  of  this  work  are  presented  here.  More 
detailed  descriptions  are  included  ir  the  attached  preprints  and  reprints. 

Ultrafast  operation  of  RPG-VCSELs: 

High  speed  pulsed  opertion  of  surface-emitting  lasers  is  necessary  for 
applications  involving  information  processing.  We  have  investigated 
experimentally  and  theoretically  the  picosecond  dynamics  of  resonant 
periodic  gain  (RPG)  vertical  cavity  surface  emitting  lasers  (VCSEL)  optically 
pumped  directly  into  the  gain  region.  The  pulsewidths  and  pulseshapes 
obtained  from  the  cross  correlation  of  pump  and  signal  pulses,  showed  a  fast 
rise  and  a  relatively  slower  fall  time.  This  was  the  first  report  of  time- 
resolved  output  pulse  shapes  for  these  lasers.  The  shortest  pulse  width  of  ~ 
13.4  pd  esd  obtained  at  a  pump  fluence  (-  5  x  threshold)  of  0.6  nJ.  The 
dependence  of  pulseshape  and  pulse  delay  agree  well  with  a  simple  rate 
equation  gain  switching  model.  This  model  gives  a  cavity  lifetime  of  8.3  ps 
and,  thus,  a  lower  bound  on  the  cavity  mirror  reflectivites  of  RiR2>0.99.  This 
is  the  first  direct  measurement  of  the  mirror  reflectivities  and  is  in  good 
agreement  with  the  values  calculated  from  simple  Fresnel  reflectivity 
models. 

Grating  Coupling  to  Surface-Plasma  Waves: 

Work  on  this  project  was  primarily  completed  and  reported  in  the  prior 
annual  report.  Papers  based  on  this  work  have  now  been  published  and  are 
included  in  the  appendix. 

High-Speed  Si  Schottky  Photodiodes: 

Relatively  little  work  has  been  reported  on  high-speed  ultraviolet  detectors. 
As  ultrafast  UV-sources  become  more  available,  these  is  increasing  interest 
in  high-speed  detectors  for  monitoring  photochemical  and  other  high-speed 
phenomena.  In  addition,  UV  is  alway  of  interest  for  space-based  systems 
because  of  the  reduced  optics  size  associated  with  the  shorter  wavelength. 
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Application  will  require  development  of  improved  detectors  as  well  as 
sources.  We  have  developed  and  tested  (UV-heterodyne  characterization  at 
334  nm)high-speed  planar  Ni-Si-Ni  Schottky  barrier  photodetectors.  For  a 
4.5-pm  gap  geometry,  a  16-GHz  3-dB  heterod3me  response,  limited  by 
packaging  effects,  is  measured  at  a  wavelength  of  334  nm.  This  is  the 
highest  reported  speed  for  a  Si  detector  at  any  wavelength.  A  photodiode 
transport  model,  including  carrier  drift,  diffusion,  and  recombination  is  in 
good  agreement  with  the  measured  detector  frequency  response.  This  model 
gives  a  24-GHz  3-dB  point  for  the  present  detector  exclusive  of  packaging 
response  and  predict  a  86-GHz  3-dB  response  for  a  1-pm  gap  geometry. 
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Ultrafast  Operation  of  Optically  Pumped 
Resonant  Periodic  Gain  GaAs  Surface 
Emitting  Lasers 

A.  MUKHERJEE,  M.  MAHBOBZADEH,  C.  F.  SCHAUS,  and  S.  R.  J.  BRUECK 


/(ftnracf— Gcncratioa  and  charactcrizatiaa  of  pkoaecoad  optical 
pabca  from  vertical<avil]r  resonani-pcriodic-taia  GaAs-AIGaAa  tmr- 
facc^mlttiag  lascn  opUcalljr  pumped  bp  picosecond  dye-iascr  pulses  Is 
reported.  The  output  pnlseshapc  was  obtained  from  tbe  cross  correia- 
Uon  of  pump  and  signal  pulses.  Dependence  of  signal  pulsewldtb  and 
pulse  delay  on  pump  power  were  invesUgated.  Tbe  results  are  in  good 
agreement  with  a  simple  rale  eqnatioB  model  of  tbe  pulse  formatloa.  A 
cavity  lifetime  of  S.3  ps,  compared  with  a  gain  medium  transil  time  of 
-0.1  ps,  is  determined  for  Ibis  very  high-(?  structure. 

I.  Introduction 

INTEREST  in  vertical-cavity  surface-emitting  lasers 
(VeSEL’s)  has  grown  recent^.  The  rapid  progress  has 
been  due  largely  to  the  development  of  high-reflectivity  Bragg 
reflectors  which  are  grown  (by  either  MOCVD  or  MBE)  as 
an  integral  part  of  the  laser  structure.  ITie  relative  simplicity 
of  fabrication,  the  potential  for  low  threshold,  high  efficiency 
lasers,  and  for  high  density  two-dimensional  array  operation 
are  especially  attractive,  ^veral  alternate  configurations  of 
the  gain  region  including  simple  beterostructures  [1],  multi¬ 
ple  quantum  weUs  [2],  resonant  periodic  gain  (RPG)  multiple 
quantum  wells  (3],  (X/2  separation),  and  single  quantum 
well  [4],  have  been  demonstrated  with  either  elecdical  or 
optical  excitation. 

For  applications  involving  information  processing,  high¬ 
speed  pulsed  operation  of  these  lasers  is  necessary.  The  very 
short  cavi^  length  (-2-5  ftm)  results  in  a  very  rapid  photon 
transit  time;  however,  the  very  high-Q  resonators  necessary 
for  laser  operation  substantially  lengthen  the  cavity  lifetime, 
and  hence  the  gain-switched  pulsewidth.  Optical  pumping 
with  ultrafast  laser  pulses  provi^  an  important  technique  for 
evaluating  the  fimdamentid  laser  temporal  behavior  without 
die  conqilications  of  circuit  panmtics  which  inevitably  effect 
electrical  pumping  measurements.  A  modulation  bandwidth 
of  8  GHz  has  been  repotted  by  Jewell  e/  al.  [5].  Raja  et  el. 
[6]  reported  optical  pumping  results  for  a  2(>-period  RPG 
GaAs-GaAlAs  laser  excited  at  620  nm  where  the  dominant 
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excitation  was  into  the  first  Bragg  reflector  followed  by 
photoluminescence  pumping  of  the  gain  region.  Time  delays 
between  the  pump  pulse  and  the  RPG  output  were  ~  1  ns 
dominated  by  the  photoluminescence  lifetime;  instrumenta¬ 
tion-limited  pulse  width  measurements  were  - 1  ns.  Re¬ 
cently,  Karin  et  al.  [7]  reported  on  autocorrelation  measure¬ 
ments  for  a  optically-pumped  MBE-grown  20-period  RPG 
VeSEL  with  30-nm  quantum  wells.  Pulsewidths  of  <  6  ps 
were  reported  for  a  comparatively  low-Q  cavity  (92.2  and 
98%  mirror  reflectivities). 

In  this  letter,  we  report  on  experimental  observation  and 
theoretical  simulation  of  the  picosecond  dynamics  of  RPG 
VeSEL’s  optically  pumped  directly  into  the  gain  region.  Tbe 
pulsewidths  and  pulseshapes  obtained  from  the  cross  correla¬ 
tion  of  pump  aixl  signal  pulses,  showed  a  fast  rise  and  a 
relatively  slower  fall  times.  This  is  the  first  report  of  time-re- . 
solved  output  pulse  shapes  for  these  lasers.  The  shortest^ 
pulsewidth  of  - 13.4  ps  was  obtained  at  a  pump  fluence 
(approximately  five  times  threshold)  of  -  0.6  nJ  (4-ps  pump 
pulsewidth  in  a  -  IO-/1  diameter  spot).  This  same  structure 
was  previously  used  for  CW  investigations  [3]  and  had  very 
high  calculate  mirror  reflectivities  (R1R2  -  99.6%).  The 
dependence  of  the  output  pulsewidth  and  delay  on  the  pump 
power  are  reported.  The  pulsewidth,  pulseshape,  and  delay 
agree  well  with  a  simple  rate  equation  model. 

n.  Theory 

The  normalized  rate  equations  describing  die  temporal 
behavior  of  the  carrier  density  and  the  optical  intensity  in  the 
laser  cavity  are  given  by 

(1) 

dS 

—  =  -S  +  S(n-/,J  (2) 

where 

r  =  time  normalized  to  the  cavity  lifetime  (=  f/r,) 
n  =  carrier  density  (normalized  to  the  clamped  carrier  den¬ 
sity  above  threshold,  neglecting  transparency,  n,  = 

1  /ffr,  where  is  the  gain  coefficient  per  carrier) 

5  =  photon  density  (also  normalized  to  n^) 

P  =  optical  pumping  rate  (normalized  to  the  cw  threshold,  g 
neglecting  transparency  requirements,  P^  =  1//St^t,)  V, 
T,  =  carrier  lifetime  in  the  gain  region 
T,  =  cavity  lifetime,  and 
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Fig.  1.  Time  evolution  of  pump  pulse,  carrier  populatioD  and  photon 

output  based  on  a  simple  rate  equation  model  of  a  gain-switched  semicon¬ 
ductor  laser. 

n„  =  the  carrier  concentration  (normalized  to  required 
for  transparency. 

For  a  pump  pulse  with  a  Gaussian  temporal  profile  of 
FWHM  -  4  ps,  (1)  and  (2)  were  numerically  integrated 
using  a  simple  predictor-corrector  algorithm.  The  result  for 
the  temporal  evolution  of  the  carrier  and  photon  density 
ccHTesponding  to  a  pump  power  of  5  is  shown  in  Fig.  1. 
The  time  scale  is  normalized  to  the  cavity  lifetime  and  the 
vertical  axis  to  n,.  The  pump  pulse  is  also  shown  in  this 
figure.  Notice  that  the  carrier  density  rises  as  the  integral  of 
the  pump  pulse  and  then  undergoes  a  very  slow  spontaneous- 
emission  decay  while  the  [rfioton  density  is  building  up  frmn 
very  low  initial  levels.  The  gain  switched  output  pulse,  along 
with  a  concomitant  rapid  depletion  of  the  carrier  density 
occurs  after  a  delay  of  several  cavity  lifetimes.  Note  that  the 
output  pulse  is  asymmetrical  with  a  relatively  fast  rise  and  a 
slower  decay.  Following  the  output  pulse,  the  carrier  concen¬ 
tration  is  depleted  to  below  iimhold  and  decays  with  the 
spontaneous  lifetime.  This  model,  neglecting  gain  and  ab¬ 
sorption  saturation.  Auger  recombination,  etc.  is  a  simplifi¬ 
cation  of  the  actual  situation.  However,  it  provides  a  very 
good  fit  to  the  experimental  dau.  A  simple  model  extension 
to  empirically  include  gain  saturation  has  been  presented  by 
Morton  et  al.  (9]. 

m.  Experiment 

The  VCSEL  structure  was  grown  by  MOCVD  on  a  GaAs 
substrate.  The  laser  structure  consisted  of  a  30j-period 
Gao75Alo.2sAs~AlAs  X/4-reflector  stack  (calculated  reflec¬ 
tivity  99.9%),  a  20-period  RPG  region  with  10-nm  GaAs 
quantum  wells  and  Gao,Alo.2As  X/2-spacers,  and  a  20- 
period  Gao.7sAloj3As-AlAs  X/4-ouqiut  reflector  (calculated 
reflectivity  99.7%).  The  A)  concentration  in  the  gain  region 
spacers  was  lower  than  that  in  the  mirror  stacks  to  allow  for 
direct  optical  pumping  of  the  gain  region.  Details  of  the 
MOCVD  grov^  and  laser  structure  have  been  reported 
previously  [3]. 

The  VCSEL  was  optically  pumped  with  the  ouqmt  of  a 
CW-modelocked  dye  laser  (4-ps  pulses,  82-MHz  repetition 
rate  at  740  nm).  A  typical  second-harmonic  intensity  autocor¬ 
relation  of  the  pump  pulses  showing  an  actual  pulsewidth 
—  4  ps  is  shown  at  the  inset  in  Fig.  2.  Notice  the  absence  of 
any  coherence  spike  in  the  autocorrelation  indicating  stable 
dye-laser  operation  [8].  This  pulse  train  was  focused  to  a  spot 
diameter  of  - 10  nm  using  a  10  x  -microscope  dijective 
(available  average  power  at  the  gain  region  was  60  mW 
accounting  for  optic  loss,  measured  reflectivity  from  the 
VCSEL  front  surface  and  estimated  absorption  in  the  top 


Fig.  2.  Crots-comlatioii  between  the  RPG  VCSEL  output  pulse  ind  the 
dye  laser  pump  pulse  for  a  pumpb^  level  ^iproximatefy  five  times 
threshold.  Note  the  asymmetry  of  the  output  pulx  (solid  line)  in  good 
agreemern  with  the  calculated  pulse  shape  (dotted  line).  The  theoretical 
time  delay  as  well  as  the  output  level  were  arbitrarily  adjusted  for  this  fit. 
The  insert  shows  the  auto  correlation  of  the  pump  pulse  of  -4  ps 
FWHM. 

mirror  (15%);  the  VCSEL  output  was  collected  by  the  same 
objective  and  directed  into  the  correlator  along  with  a  split-off 
portion  of  the  pump-pulse  train.  The  harmonic  crystal  was 
rotated  to  provide  a  cross-correlation  response  at  the  398-nro 
sum  frequency  of  the  740'nm  pump  and  the  862-nm  VCSEL 
output.  This  procedure  provides  the  pulse  width  and  a  rela¬ 
tive  pulse  delay.  For  the  shortest  pulse  delay,  both  the  pump 
pulse  reflected  from  the  VCSEL  and  the  output  pulse  were 
observed  with  a  fast  detector  and  oscilloscope  combination  to 
provide  a  calibration  of  the  peak-to-peak  time  delay  ( - 16 
GHz  bandwidth  to  determine  the  absolute  delay);  this  intro¬ 
duces  a  relative  error  of  ~  ±  3  ps. 

The  cross-correlation  signal  for  the  shortest  VCSEL  output 
observed  is  shown  in  Fig.  2.  The  pulse  width  is  -  14  ps  with 
a  ^  rise-  and  a  relatively  slower  fall-time.  Taking  the  finite 
pump  pulsewidth  into  account,  the  pulsewidth  of  the  VCSEL 
was  estimated  to  be  -  13.4  ps.  The  small  prepulse  is  due  to 
a  residual  satellite  pulse  present  in  the  pump  pulse  train  used 
for  this  particular  measurement.  The  r^tively  longer 
pulsewidth  from  our  VCSEL  compared  with  a  FWHM  <  6 
ps  reported  in  [7]  is  attributed  to  our  higher  mirror  reflectivi¬ 
ties,  resulting  in  longer  cavity  lifetimes.  The  dotted  line  is  the 
numerical  result  from  the  rate  equations  for  the  parameten, 
Tp  =  8.3  ps,  =  0.5,  T,  =  200  T,,  and  P  =  5.  The  time 
delay  and  pulse  amplitude  were  independently  adjusted  to 
match  the  experimental  result.  Good  qualitative  agreement 
with  the  details  of  the  experimental  pulsesh^  is  obtained. 

An  input-ouqxit  curve  for  the  ultrafast  operation.  Fig.  3, 
of  this  VCSEL  showed  a  threshold  of  10-mW  average  power. 
The  external  slope  efficiency  of  the  theoretical  input-output 
curve  was  adjusted  to  the  experimental  value  of  ~  10%.  The 
cross-correlation  data  showing  the  pulsewidth  and  peak-to- 
peak  delay  as  a  function  of  pump  power  is  shown  in  Fig.  3. 
The  x-axis  is  pump  power  in  units  of  threshold  power  and 
the  y-axis  is  time  in  units  of  cavity  lifetime.  The  solid  lines 
represem  the  theoretical  fit  for  the  identical  parameters  used 
to  fit  the  pulseshiqie  (Fig.  2).  In  fitting  these  dau,  the 
calculated  results  for  P  were  renormalized  by  1  /(I  +  n„)  to 
account  for  die  transparency  requirements.  That  is,  the  defi¬ 
nition  of  neglecu  transparency  requirements;  since  the 
experimental  daU  are  normalized  to  the  actual  threshold, 
including  transparency  requirements,  the  model  values  of  P 
are  divided  by  (1  +  n„)  for  comparison.  Again,  there  is  very- 
good  agreement.  There  exisu  some  uncertainty  in  the  r,  and 
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Fir.  3.  P«>wer  dependence  of  the  RPG  VCSEL  output  polsewidth,  delay 
from  the  pun^  pulae,  and  output  power.  The  power  axis  haa  been 
normalized  to  the  observed  threshold  power  for  the  experimental  results 
and  to  the  cw  threshold  (including  transparency  rerpiireinents)  for  die 
model  resultt.  Good  agreement  it  obtained  for  a  cavity  lifetime  of  8.3  ps. 
The  dashed  curve  associated  with  the  pulse  delay  data  is  an  approximate 
analytic  result  (see  text). 

values;  comparable  fits  were  obtained  for  a  relatively 
wide  range  of  values.  A  time-resolved  photoluminescence 
experiment  to  independently  measure  r,  U  planned;  this  will 
allow  a  more  precise  estimate  of  from  this  model.  The 
experimental  points  for  both  pulsewidth  and  delay  deviate 
firom  the  theoretical  fit  towards  longer  times  for  higher  pump 
powers.  This  may  result  from  gain  saturation  or  spectral  hole 
burning  leading  to  absorption  saturation  of  the  pump.  Any 
mechanism  that  decreases  the  gain  below  the  linear  regime  of 
the  model  will  result  in  a  saturation  in  die  pulsewidth  and 
delay  above  that  predicted  by  this  simple  model. 

IV.  Discussion 

It  is  informative  to  compare  the  cavity  lifetime  with  a 
simple  cavity  transit  time.  The  cavity  consists  of  the  20  X/2 
periods  of  tte  gain  region  plus  the  distributed  mirror  layers. 
An  effective  cavity  length  is  given  by  including  the  mirror 
1/e  lengths  (~  10  X/4  sections  for  each  mirror)  along  the 
gain  regkm;  this  gives  round-trip  transit  time  r,  of  ~  9  x 
10~  s.  The  effective  cavity  parameters  are  related  to  these 
times  by  the  simple  relationship 

1  -  =  t,/t,  (3) 

where  X,  and  are  the  mirror  reflectivities  and  a  is  the 
(effective)  loss  coefficient  due  to  background  and  mirror 
fcaHvring  and  free  carrier  absorption.  Assuming  that  the 
mirror  reflectivities  are  the  theor^kal  values  of  99.9%  and 
99.7%  for  the  present  structure,  (3)  can  be  used  to  estimate 
_  0.993.  This  result  is  consislent  with  the  observed 
differential  quantum  efficiency  in  this  eiqieriment  of  ~  12%, 
but  is  inconsistent  with  the  much  higher  differential  (and 
indeed  overall)  power  conversimi  efficiencies  observed  for 
CW  excitation  in  this  same  structure  [3].  This  discrepancy 
may  be  due  to  higher  photon  losses  in  the  pulse  exciution 
case  where  much  higher  intensities  leading  to  saturation  and 
muhiphoton  excitation  processes  are  possible.  Detailed  mea¬ 
surements  as  a  function  of  mirror  reflectivity  are  underway  to 
rest^e  this  issue. 

It  is  possible  to  obtain  an  approximate  analytic  result  for 
the  time  delay  between  the  exciution  and  output  pulses, 
itamiming  that  both  pulses  are  short  compared  with  the  time 
delay  from  (1)  and  (2).  During  this  time  delay,  the  photon 
density  is  negligible  so  that  (1)  simply  gives  an  exponential 
decay  for  the  carrier  density  from  du  initial  exciution  pulse 


at  the  carrier  lifetime  r,.  This  result  is  then  substituted  back 
into  (2)  which  can  be  integrated  to  yield 

S(t)  *  SoeM»{-T(l  +  nj  +  £(t,/t,)(1  -  e""r/Tj} 

(4) 

for  a  delay  tune  much  less  than  r^.  Here,  £  is  the  integrated 
energy  in  the  pump  pulse.  The  second  exponential  in  (4)  can 
be  expanded  with  the  result  that  for  an  exponential  growth  of 
A  for  lasing  the  very  simple  result 

r^-Arj[P-{\+n^)\  (5) 

is  obtained.  TraditionaUy,  values  of  A  between  20  and  30 
have  been  invoked  for  amplification  above  the  background 
levels,  a  value  of  23.S  provides  the  best  fit  with  the  detailed 
numerical  calculation.  This  result  is  plotted  as  the  dashed 
curve  in  Fig.  3.  Good  agreement  is  obtained  with  the  more 
exact  numerical  calculation. 

We  have  reported  generation  and  characterization  of  ultra¬ 
fast  optical  pulses  from  a  GaAs-GaAlAs  RPG  VCSEL  opti- 
'  ally  pumped  by  4-ps  duration  dye-laser  pulses.  Cross-corre¬ 
lation  experiments  yield  a  minimum  pulsewidth  of  ~  13.4  ps 
as  well  as  the  asymmetric  pulse  shape,  and  the  power  depen¬ 
dence  of  the  pulse  duration  and  delay.  These  results  are  in 
good  agreement  with  a  simple  rate  equation  model  of  the 
gain-switched  pulses.  Analysis  of  the  rrsults  point  to  varia¬ 
tions  in  cavity  losses  between  CW  and  pulsed  operation  of 
these  lasers.  Clearly,  a  shorter  cavity  length  and/or  lower 
minor  reflectivities  will  yield  shorter  pulses.  Further  investi¬ 
gations  are  underway. 
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A  Simple  High-Speed  Si  Schottigr  Photodiode 

B.  W.  MuUinf,  S.  F.  Soares,  K.  A.  McArdle,  C.  M.  ^KHIsoq,  aasd  S.  R.  J.  Biueck 


yf  ftunKf— DMlga,  Mbricalioa,  tad  UV-hcterodyne  chsncter- 
Isatloa  of  N1-8I-NI  ■ctal-tenkoadactor-aMtal  Schottky  bar- 
rtor  pkotodctMton  Is  nportod.  Ptaaar  detectors  were  febrl- 
CRled,  wHk  •  riBMde  3-level  ■thography  proctas  oa  balk  9,  la 
both  teple-gap  aad  laterdigitated  gcoaietiiea  wHk  gap  diaiea- 
iloas  firoai  1  to  5  pa.  IVeqaeaey  reapooa  of  thea  devloa  was 
characterlad  aaiag  a  CW-hner  heterodyac  tystca  at  334 J  aa. 
For  a  4  J-pa  laterAgItated  devict,  a  3-0  respoaa  of  Id  GHt  is 
aeasared,  gMag  23  GHx  wbea  decoavohred  froa  the 
pacfcagc/coaaeetor.  A  detailed  theoretical  aodd  of  the  photodi¬ 
ode  respoaa  tacorporatlag  effects  of  cankr  traasport  aod  de- 
via  geoasetry  is  la  czccDeat  agreeaaat  with  Oe  oicasareaeat 
This  aodal  predicts  aa  S6-GHa  3-dB  respoaa  fa  the  1-pa  gap 
geoasetry  drideea. 


iKTHODUCnON 

HERE  has  been  exlatsive  research  on  high-speed  low- 
ooise  pbotodetectors  in  the  0.8-1.5  pm  wavelength  range 
for  use  in  optical  communication  systems.  Four  general 
device  designs,  p-i-n  [1],  (2],  vertical  Schottky  barrier  [31, 
planar  Schottky  barrier  [4]-[7],  aod  avalanche  fdiotodiode 

[8] ,  have  been  developed.  Device  materials  include  Si  [S], 

[9] ,  InGaAs  [8],  and  strained-layer  superlattices  [10].  In  the 
near  infrared  regkm,  v«y  fast,  highly  responsive  detecton 
have  been  fabricated.  Frequency  refuses  up  to  67  GHz  and 
gpiantam  efficiencies  of  80%  Ittve  been  reported  for  p-i-n 
devices  [2].  The  ffistest  Sebottky-barrier  devices  r^orted 
have  a  llO^Hz  frequency  response  [4]. 

In  contrast,  relatively  little  work  has  been  reported  on 
ultraviolet  detecton.  As  ultrafast  UV  sources  become  avail¬ 
able,  diere  is  increasing  interest  in  high-qwed  detectms  fw 
monitoring  photochemical  and  odier  high-qieed  phenomena. 

The  most  salient  materials  duuacteristic  at  UV  wave¬ 
lengths  that  determines  detector  design  is  the  very  high 
absorption  coefficient  of  ~  10*  cm  for  almost  all  semicon- 
dudon.  This  dictates  a  design  with  very  shallow  junctions  in 
traditional  vertical  geometries  and  fiivon  Ae  (danar  struc¬ 
ture. 
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High-speed  UV  detectm  hive  been  rqiorted  using  vertical 
Schottky  barrier  [11],  planar  interdigitated  Schottky  barrier 
[9],  Si  p-n  junction,  and  vertical  p-i-n  with  a  thin  p-layer 
[12],  wiA  ody  moderate  success.  A  promising  design  used  a 
pla^  interdigitated  School^  barrier  widi  a  low-frequency 
quantum  efficiency  of  12%  at  337  nm  and  a  fie^ncy 
response  of  5.3  GHz,  measured  at  800  nm  [9).  No  picosec¬ 
ond  impulse,  autocorrelation,  or  heterodyne  teq>onse  meas¬ 
urements  at  UV  wavelengths  have  been  rqwited. 

We  report  on  fsbricatitm  and  UV-heterodyne  characteriza¬ 
tion  of  high-speed  planar  Ni-Si-Ni  Schottky  barrier  pbotode- 
tectors.  For  a  4.5-/tm  gap  geometry,  a  16-GHz  3-dB  hetero¬ 
dyne  response  limited  by  packaging  effects,  is  measured  at  a 
wavelengtii  of  334  nm.  This  is  the  highest  reported  speed  for 
an  Si-based  detector.  A  photodiode  transp(»t  model,  includ¬ 
ing  carrier  drift,  diffusitm,  and  reemnbination,  is  in  good 
agreement  wiffi  the  measured  detector  frequency  re^xMise. 
This  model  gives  a  24-GHz  3-dB  point  for  Ae  i»eseat 
detector  exclusive  of  packaging  limitations.  For  a  l-/im  gqi 
geometry,  the  model  predicts  an  86-GHz  3-dB  response  for 
these  sinqile  Si  detectors. 

The  detector  consists  of  a  side^^-side  pair  of  metal- 
semicondoctor  junctions  comprising  a  metal-semicon- 
ductw-metal  (MSM)  Aode.  Ibe  exposed  semiconductor, 
between  the  metal  pads,  is  folly  dqileted  wiA  an  appre^niate 
applied  bias.  This  dqileted  regkm  is  die  photosensitive  region 
of  die  Aode,  and  hi^  applied  Sdds  result  in  very  fast, 
saturated  drift  velocities  of  tte  charge  carriers,  ft  is  important 
to  ^ipreciate  that  boA  electron  and  htde  velocities  are  equally 
irnportam  in  this  configuiatioo.  In  this  regime.  Si  offers 
perftHmance  comparable  to  GaAs  and  InP,  in  contrast  to  the 
situation  in  which  doly  electron  velocities  and  mobilities  are 
impmtant.  The  direct  access  to  the  depletion  region  by  UV 
photmis  avoids  the  short  absorption  lengA  problem  inherent 
in  other  designs. 

A  common  pbnomenon,  rqxitted  fK  all  MSM  devices,  is 
a  tporious  gain  meefaaniam  kaffing  to  an  anomalous  effective 
quantum  efficiency,  >  100%,  at  low  frequencies  [9],  [13]. 
Tentative  explanations  have  been  pot  forth  [14],  but  no 
definitive  models  have  been  (HPesented.  This  {Aenomenon  is 
also  seen  in  the  devices  terted  in  diis  study.  Our  results 
indicate  that  this  gain  is  most  likely  photoconductive  resulting 
from  junction  leak^;  a  mme  detailed  analysis  will  be 
publisM  elsewhoe. 

Fabrication  of  Ni-Si-Ni  MSM  Soiama  Barriex 

PHOTtHHODBS 

MSM  jAotodiodes  were  fiArkated  wiA  Ni  Sdiottky  con¬ 
tacts  on  bulk  Si  (n  —  8  X  10'*  cm~’).  A  pbuur,  latnal 
geometry  wiA  boA  ccmtacts  on  die  top  surfooe  of  die  Si  was 
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Flf.  I.  MSM  dettctor  ftraclure.  Boodiof  pid:  75  x  75  Wiadow: 
25  X  20  juB^  (•)  Cr/An  (200  nm}-  (b)  N1  (30  nm).  (c)  SiOi  (W  obi),  (d)  SI 
(bulk). 

used;  both  simple-gsp  and  interdigitated,  or  seipentioe-gi^, 
devices  were  fabricated  with  g^>  widths  varying  from  1  to 
4.S  iim.  A  representative  structure  is  shown  in  Fig.  1.  The 
first  layer  is  an  80-nm*thick  SiO]  film.  Windows,  ranging 
from  100  to  SOO  ittB?,  expose  the  semiconductor  surface 
where  the  Schottky  contacts  are  formed.  The  oxide  serves  to 
iaolMf  optically  inactive  metalluation  areas  and  reduce  June* 
tion  capacitance  and  dark  current.  The  second  layer  is  a 
30-nm-thick  Ni  film  with  a  series  resistance  of  SO-M  Q  and 
a  transmission  of  17%  at  334.S  nm.  The  bTi-Si  contacts  form 
the  MSM  photodiode.  The  photosensitive  semiconductor  re¬ 
gion  between  die  contacts  varies  in  area  torn  IS  to  90  /im^ 
in  the  simple-gap  geometry.  The  final  layer  is  a  30/170  nm 
Cr/Au  film  to  form  the  75  x  7S  /im*  bonding  pads. 

Photodiodes  were  fabricated  using  cmiventional 
toUthography,  including  standard  liftoff  techniques  for  metal 
definition.  Si  wafers  were  cleaned  and  insert^  into  a  ptt- 
heated  furnace  (1 100*Q,  with  an  O]  flow  of  4 1/min,  for  2S 
min,  to  grow  the  80-nm  oxide.  The  windows  in  the  oxide 
film  were  etched,  and  the  metal  films  dqx»ited  in  three 
photolithographic  steps.  Metal  films  were  deposited  in  an 
electron-beam  evqx>rator  at  a  rate  of  O.S  nm/s. 

After  fabrication,  wafers  were  cleaved  to  separate  dies 
irKorporating  detectors  with  variations  in  geometry  and  di¬ 
mensions.  An  individual  die  was  mounted  at  the  edge  of  an 
Au-on-alumina  microstrip  high-speed  package  using  heat- 
conductive  epoxy.  The  package  was  mounted  on  an  Au-coated 
A1  block  arxl  an  SMA  connector  (18-GHz  barxlwidth)  was 
pressure  fitted  to  the  strip.  Gold  wires  (2S-fim  diameter)  we.'e 
wire  bonded  between  the  photodiode,  die  package,  and 
ground.  Care  was  taken  to  ensure  short  bond  wires  so  that 
pariraging  inductance  did  not  affect  die  measuremeats.  Darit 
currents,  in  this  nonoptimized  design,  were  <50  nA  for  a 
25-;tm-long  device  at  a  10-V  bias.  For  broad-band  photodi¬ 
ode  applications,  the  noise  associated  with  this  low  dark 
curirat  is  negligible  in  comparison  to  RF-amplifier  input 
nmae. 

EtCPEMMEKTAL  AMtANaa«fENT 

A  UV-heterodyne  systmn  was  developed  to  characterize 
these  detecton  115].  A  single-mode  Ar*-ioo  laser,  X  =  334.5 
ran,  and  an  intracavity  (kxiUed  (LflO,),  frequency  stabilized 
dye  laser  provided  the  local  oscillator  and  signal  beams, 
le^ectively.  Careftil  attettkm  was  placed  on  beam  and  polar¬ 
ization  overlap.  A  typical  local-oscilbtor  intensity  tff  10'’ 
W/cm’  (-  2-/am  FWHM)  tf  the  detector  was  necessary  to 
re^  shot-otdse-limited  detection.  The  gain  of  the  RF- 


Pif.  2.  Picqueney  itqxxue  of  a  4.5-Mm  gap  intetdigitaied  ddector  at  30  V 

bias. 

anqilifiers  and  cabling  fitun  the  SMA  connector  to  the  RF- 
spectrum  analyzer  was  calibrated  from  1  MHz  to  18  GHz 
using  a  vector  network  analyzer.  Absolute  lesponsivity  R(.a) 
was  confirmed  by  conqiarison  to  the  measured  dc  responsiv- 
ity  of  a  well-characterized  low-frequency  commercial  detec¬ 
tor. 

Results 

Measuranents  of  /ffoi)  were  made  on  several  different 
detectors  under  a  variety  of  conditions.  The  bulk  of  these 
measurements  were  used  to  validate  an  MSM  detector  re¬ 
sponse  nxidel  that  predicts  high-speed  detector  performance 
as  a  function  of  bias  and  beam  diameter  and  position  on  the 
detector,  using  Si  carrier  transport  prtqierties  [15].  The  model 
includes  effects  of  carrier  drift,  dif^ion,  and  recombination. 
A  mote  detailed  description  will  be  published  elsewhere.  The 
model  was  compered  to  the  measured  frequeiKy  response  for 
g^)  dimensirms  of  4.5,  3,.and  1  ^m  with  differing  biases  aixi 
bew  positions. 

curves  are  shown  for  a  4.5-;im  gap  interdigitated 
detector  at  30-  and  lO-V  bias.  (Figs.  2  and  3)  Heterodyne 
signal  levels  were  greater  than  20  dB  above  the  shot-noise 
floor  across  the  entire  frequeix^  range;  the  —  ±1  dB  varia¬ 
tions  are  doe  to  package  and  electronics  resonances  as  can  be 
seen  by  noting  the  pattern  similarity  in  the  two  figures.  The 
l-fim  FWHM  beam  was  centered  in  the  detector  gap,  and 
R(,u)  was  measured  out  to  18  GHz.  The  30-V  data  (Fig.  2) 
have  been  extended  to  show  die  low-frequency  gain  r^on, 
inherent  in  all  MSM  photodiodes.  Two  theoretical  plots  are 
shown.  The  upper  curve  in  both  grafAis  is  the  detector  model 
which  is  in  good  agreetneot  witii  die  data  out  to  about  14 
GHz.  At  these  frequeodes,  the  bandwiddi  limitations  of  the 
SMA  connector  on  the  detector  package  become  importam. 
The  SMA  connector-detector  package  can  be  mocfeled.as  the 
product  of  a  sirrqile  single  p^,  with  a  characteristic  fre¬ 
quency  of  18  G&,  due  to  die  SMA  connector,  and  die 
detector  reqxxise.  This  leads  to  the  lower  modeling  curve. 
An  excellent  fit  is  obtained  out  to  the  eiqierimental  frequency 
limit  of  18  GHz,  constrained  by  available  amplifiers.  The 
effect  of  die  SMA  connector  can  be  seen  on  both  the  30-  and 
lO-V  data.  The  variation  in  the  frequency  response  for  these 
two  cases  arises  because  of  the  (slight)  change  in  carrier 
velocities  at  die  two  applied  field  strengths  in  this  saturated 
regime.  The  ody  adjutfaUe  parameter  is  die  overall  quantum 
efficiency. 

An  extrapolation  of  the  detector  bandwidth  can  be  made 
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Fig.  3.  Frequency  response  of  a  4.5-Mm  gap  interdigiiated  detector  at  10  V 

bias. 


Rg.  4.  Theoretica]  frequency  response  curves  for  varying  gap  dimeiisioos. 
The  calculation  was  carried  out  for  a  fixed  2-ttm  iHuminatioa  spot  and  a 
fixed  biu  of  10*  V/cm.  The  reduced  responsivity  for  die  smaller  detecton 
is  simply  a  geometrical  eflect  due  to  the  reduced  overlap  between  the  beam 
and  the  active  detector  area. 

from  the  model,  assumii^  an  ideal  RF  package.  At  30-V 
bias,  the  bandwidth  of  diis  4.S-/tm  gap  det^or  is  21 .5  GHz. 
Increasing  the  bias  wUl  result  in  slighdy  higher  carrier  veloci¬ 
ties  and  should  give  a  24-GHz  bandwidth  for  a  4.S-^m  gap 
Si-based  detector.  Scaling  the  response  to  smaller  detector 
gap  widths  of  3  and  1  iim,  the  moM  predicts  bandwidths  of 
34  and  86  GHz,  respectively  (cf.  Fig.  4).  The  lower  respon¬ 
sivity  shown  for  the  smaller  detectors  is  simply  a  geometrical 
effect  as  the  calculation  is  carried  out  for  a  incident  spot 
size.  This  high-speed  response  is  con^Miable  to  the  fastest 
rqKHted  GaAs  detectors.  The  high-fre^ncy  ultraviolet  re- 
qxHisivity  of  these  detectors,  0.032  A/W  (—  12%  quantum 
efficiency)  with  no  antireflective  coating,  is  very  high  for 
detectors  with  diese  large  bandwkfeha.  The  quantum  effi¬ 
ciency  of  the  smaller  gap  detectors  scales  only  with  die  gap 
dimension/spot  sue  ratio.  Thus,  if  the  incident  spot  is  smaller 
than  the  gap  dimension,  then  the  quantum  efficiency  remains 
nearly  constant.  Smaller  detectors  and  tighter  focusing  will 
lead  to  a  significandy  iffl{HOved  quantum-efficiency  band¬ 
width  product. 

Conclusion 

A  high-speed  MSM  photodiode,  based  on  bulk  Si  and 
fobricated  by  a  sinqile,  inexpensive  process,  has  been  denxni- 
strated.  The  bandwiddis  of  these  devices  compare  well  to  die 
fostest  reported  GaAs  devices.  A  major  advanuge  of  this 
device  for  UV  apfdications  is  the  relative  robustness  of  Si 
oonqMred  to  GaAs  when  subjec  d  to  intense  ultravidet 
radiation.  GaAs  devices  were  also  fobricated;  however,  they 
showed  severe  degradation,  on  a  time  scale  of  minutes,  when 
subject  to  the  intense  local  oscillator  irradiation.  This  degra¬ 


dation  was  irreversible  and  severely  limits  dte  apfdkability  of 
the  unpassivated  GaAs  devices  for  ultraviolet  sp^cations.  In 
cmitrast,  the  Ni-Si  devices  were  extremely  robust  widi  no 
changes  in  speed  or  responsivity  observed  over  time  scales  of 
several  months— both  for  devices  stored  in  laboratory  ambi¬ 
ent  and  devices  subjected  to  many  hours  of  high-intensity  UV 
irradiation. 

The  simplicity,  manufacturability,  interconnection  compat¬ 
ibility,  and  robustness  of  the  planar-Si  MSM  fdiotodiode 
nudce  it  a  strong  candidate  for  plications  where  t^  ultimate 
connection  is  to  a  silicon  VLSI  circuit.  The  high-speed 
high-iesponsivity  detector  presented  here  is  an  ideal  candi¬ 
date  for  many  plications  due  to  its  materials  compatibility 
with  standard  Si  technology  and  the  absence  of  (xnnplex  and 
epnsive  q>itaxial  growth. 
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A  aystematic  experimental  and  theoretical  stud7  of  first-order  grating  coupling  to  surface  plasma  waves  exist¬ 
ing  at  an  air-Ag  interface  ia  presented.  The  experiment  extei^  beyoird  previous  work  to  grating  depths  com¬ 
parable  with  the  grating  period.  The  grating  profiles  range  from  sinusoidal  to  rectangular.  For  TM-poIarized 
incident  radiation  this  grating  depth  range  includes  the  entire  spectrum  of  surface-plasma-wave-radiation 
coupling— from  underdan^ied,  to  nearly  100%  cou{ding,  to  overdamped  and  the  disappearance  of  the  resonance 
from  Um  zero-order  reflectance  measurements.  Strong  polarization  and  absorption  effects  are  observed  for 
the  deepest  gratings.  A  simide  theoretical  model,  based  on  the  Raylei^  hypothesis  and  retaining  only  reso¬ 
nant  diffraction  terms  without  a  small-signal  approximation  being  made,  provides  good  agreement  with  the  ex¬ 
perimental  results. 


INTRODUCTION 

The  stu^y  of  the  interaction  of  light  with  periodic  struc¬ 
tures  (gratings)  on  metals  has  a  long  and  distinguished 
history.  In  1902  Wood*  first  noted  the  anomalous  behav¬ 
ior  (christened  IK^rod’s  anomalies)  displayed  diffraction 
gratings  of  large  and  rapid  changes  in  diffraction  intensi¬ 
ties  for  small  angular  and  spectral  variations.  Rayleigh*^ 
presented  in  1907  the  first  theoretical  explanation  of  these 
anomalies  in  suggesting  that  such  behavior  was  due  to  the 
cutoff  or  the  aniearance  of  a  new  spectral  order.  In  1941 
Fnno  first  distinguished  between  two  types  of  Wood 
anomalies;  (1)  an  edge  anomaly,  with  a  sharp  behavior 
related  to  the  passing  off  of  a  diffraction  order  (i.e.,  a  dif¬ 
fraction  order  passing  over  the  horizon,  9(r  to  the  surface 
normal)  and  (2)  a  resonance  anomaly  due  to  excitation  of  a 
bound  w  surface  wave  at  the  metal-dielectric  interface.* 

Surface  plasma  waves  (SPWs)  are  TM  modes  of  the 
electromagnetic  field  bound  to  the  interface  between  a 
metal  and  a  dielectric.  The  condition  for  existence  of  the 
SPW  mode  is  <  -e^,  where  tm'  is  the  real  part  of  the 
metal  dielectric  constant  and  is  the  dielectric  constant 
of  the  dielectric.  Related  modes,  first  investigated  fay 
Sommerfeld,*  exist  when  one  of  the  media  is  higUy  lossy. 
The  SFW  {diase  velocity  is  less  than  the  light  vekKity  in 
the  dielectric,  and  phase  matching  between  an  incident, 
freely  propagating  wave  and  the  SPW  is  accomplished 
eithm'  by  a  prism*  or  by  grating  coupling  techniques. 

The  SPW  di^iersion  relation  for  a  {danar  metal-air  in¬ 
terface  is  simi^  given  hy^ 

*aew  -  *o[«-/(e-  +  1)]“.  (1) 
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where  ke  -  2fl/A  is  the  free-space  optical  wave  vector. 
The  phase-matching  condition  for  excitation  of  SPWa  is 
satisfied  vdienever  *o  Bin  R  the  component  of  ko  along  the 
metal-air  interface,  satisfies  the  condition 

ko  sin  9  =  ±k8em'  +  n2iTfd,  (2) 

adiere  <f  is  the  grating  period,  9  is  the  angle  of  incidence, 
n  =  ±1,  ±2, ...  is  the  coupling  order,  and  kspm’  is  the  real 
part  of  the  SPW  wave  vector.  In  1^.  (2)  the  choice  of  a 
negative  sign  preceding  kgn'  correqwnds  to  a  SPW  mov¬ 
ing  in  the  Erection  oiq;>08ite  the  incident  wave.  This  ex¬ 
pression  assumes  that  the  grating  wave  vector  is  in  the 
plane  of  incidence;  i.e.,  the  grating  lines  are  perpendicular 
to  the  incident  wave  vectmr. 

The  SPW  dispersion  relationship  is  plotted  in  Fig.  1  for 
a  lossless,  free-electron  metal  (c.  »  1  -  «//»*),  where 
the  axes  are  nmnalized  to  the  jdasma  frequengr  w,  and 
the  corresponding  optical  wave  vector  (kf  »  w,/h).  Also 
shown  as  two  vertical  dashed  lines,  ccoresponding  to  the 
n  =  ±1  orders  in  Eq.  (2),  are  the  wave  vectors  of  a  surface 
grating.  Finally,  the  wave  vectors  accessible  in  the  n  s  o 
and  n  s  ±1  orders  by  variations  in  the  an^e  of  incidence 
are  shown  as  horizontal  lines.  This  figure  was  drawn  for 
a  grating  period  smaller  than  the  optical  wavelength 
(A/i  >  1);  note  that  there  is  only  one  point  that  satisfies 
Eq.  (2),  for  n  =  -1,  and  at  this  incident  anj^e  there  is  no 
aUowed  diffraction  order. 

In  1967  Teng  and  Stem  Hrst  detected  SPWs  optically 
by  bombarding  1200-line/mm  (833-nm-period)  gratings 
with  10-kdlf  electrons  and  observing  the  emitted  optical 
radiation.*  They  observed  that  the  line  shape  of  the 
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WAVE  VECTOR  (ky  /  kp) 

Fig.  1.  Disperaion  relation  of  SPWs  for  a  loeslesa,  free-electron 
metal  (c  1  -  «,*/«’).  The  axes  are  normalized  to  a>,  and 
k,  a  The  grating  wave  vectora  corresponding  to  n  =  ±1 

order  of  a  grating  of  period  d  are  shown  as  vertical  dashed  lines. 
The  range  of  wave  vectors  accessible  hgr  variations  in  the  irqait 
an^  from  normal  to  grazing  incidence  for  n  =  0  and  n  =  ±1  is 
shown  as  horizontal  lines.  Note  that  for  this  choice  of  parame¬ 
ters  (A/i  <  1)  there  is  (mly  one  SPW  cotipling  resonance,  and  at 
this  resonance  angle  there  is  no  propagating  diffraction  order. 

emitted  radiation  was  influenced  Iqr  the  surface  condition 
of  the  metal  but  was  independent  of  the  energy  of  the  bom¬ 
barding  electrons.  Covran  and  Arakawa  carried  out  a 
detailed  study  of  the  SPW  dispersion  curves  for  dielectric- 
metal  layers  on  concave  diffraction  gratings  and  also  de¬ 
veloped  a  quantum-mechanical  formalism  to  describe 
their  results.*  Hutley  and  co-workers  published  in  1973  a 
detailed  experimental  study  of  the  anomalies  of  sinusoidal 
profile  gratings  as  a  function  of  grating  depth^*'**  and 
characterized  the  SPW  line  shapes  for  grating  depths 
h  ranging  up  to  60  nm  for  500-nm-period  gratings  (A/d  s 
0.12).  The  integral  formalism  developed  by  Petit  and 
others'*  was  used  to  describe  these  results  with  good 
agreement  Pockrand  and  Raether,  in  an  extensive  series 
of  publkationa,'*''*  characterized  the  SPW  coiqpling  as  a 
ftinction  of  grating  period,  depth,  and  profile.  The  grat¬ 
ings  studied  were  sufficiently  deep  for  better  than  98% 
coiqding  into  the  mode  to  be  realized.  A  perturba- 
ti<m  analysis,  developed  by  Kroger  and  Kretschmann,'* 
was  used  in  modeling  these  results  with  good  agreement, 
although  clearly  the  perturbation  approach  must  break 
down  as  the  coupling  efficiency  iqiproeches  100%  A  com¬ 
plete  theoretical  treatment  of  grating  coiqiling  to  SPW’s 
was  provided  by  Mills  and  co-workers.*^  Their  apiaoach 
usee  an  integral  formulation  of  the  boundary-value  prob¬ 
lem  at  the  grating  interface  and  an  extinction  theorem 
mechanism  fidkiwing  from  Green’s  theorem  that  was  first 
iqiplied  by  Tngo  et  al.**  The  latter  stu^y  does  not  give 
simple  analytic  results  but  relies  on  extensive  compu¬ 
tational  and  numerical  evaluation.  Glass  et  al.**  and 


Weber**  simplified  this  treatment  by  developing  a  renor¬ 
malized  mode-coupling  theory  that  retained  the  resonant 
terms  and  treated  all  other  terms  within  perturbation 
theory.  Yamashita  and  Tsigi**  independently  develq;>ed  a 
differential  formulation  tiiat  treats  the  resonantly  gener¬ 
ated  fields  on  a  par  with  the  incident  fields  and  allows  for 
saturation  and  decreases  in  coupling  with  increasing  grat¬ 
ing  depths.  They  empk^ed  a  power-series  expansion  in 
kt,h  and  restricted  their  study  to  small  grating  ampli¬ 
tudes;  simple  analytic  expressions  were  obtained  for  the 
coupling  strength  as  a  function  of  h. 

In  this  pcq>er  a  comprehensive  experimental  and  theo¬ 
retical  stu()y  of  Grst-onier  grating  coupling  to  SPW’s  for  a 
wide  range  of  grating  parameters  is  presented.  The  ex¬ 
perimental  results  establish,  for  the  first  time  to  our 
knowledge,  a  relationship  between  grating  depth  and 
grating  period  for  SPW  coupling  and  extend  to  grating 
depths  that  no  longer  support  SPW’s  but  rather  show 
polarization  and  absorption  effects.  The  theoretical  iq>- 
proach  is  an  extension  of  the  Rayleigh  hypothesis  includ¬ 
ing  only  resonant  terms  in  the  Rayleigh  expansion  (cf. 
Refs.  23-25).  This  results  in  considerable  simplification; 
reasonable  agreement  between  theory  and  experiment  is 
found  out  to  depth/^iod  ratios  of  —0.5.  Gupta  et  al.^ 
used  a  similar  approach  to  describe  grating  coupling  to 
long-range  SPW  modes  on  thin,  symmetrically  bounded, 
metal  films. 

GRATING  FABRICATI(R*J 
Gratings  were  fabricated  holographically  in  positive 
toresist  layers  spun  onto  Si  substrates  with  the  use  of  the 
488-nm  line  from  a  single-mode  Ar-icm  laser.  The  details 
of  the  grating  fabrication  were  presented  elsewhere.** 
The  grating  proGles  were  approximately  sinusoidal  for 
shallow  dqiths,  evolving  towsud  rectangular  profiles  as 
the  depth  was  increased  (cf.  Fig.  4  below).  After  develop¬ 
ment  these  gratings  were  coated  with  —lOO-nm-thick 
electron-beam  evaporated  Ag  films.  Films  were  de¬ 
posited  at  room  temperature  and  at  background  pressures 
of  low  10'*  Torr. 

OPTICAL  ARRANGEMENT 
AD  of  the  measurements  reported  here  are  of  the  angular 
dependence  of  the  zero-order  reflectance  for  a  fixed- 
frequency  TM-polarized  He-Ne  laser  beam  at  633  nm. 
The  samples  were  mounted  as  one  surface  in  a  90*  comer 
refketor  attached  to  a  conqiuter-contrdled  rotation  stage. 
This  arrangement  ensured  that  the  reflected  beam  was 
always  returned  in  the  same  directiem  and  eliminated  the 
necessity  of  a  second  rotatiem  stage  for  the  detector.  The 
incident  laser  beam  was  focused  with  a  long-focal-length 
(0.5-m)  lens  to  a  spot  of  ~2  mm.  This  provided  an  angu¬ 
lar  limitation  of  0.05*,  much  smaDer  than  the  observed  re¬ 
flectance  variations.  Care  was  taken  with  the  alignment 
to  ensure  that  the  axis  of  the  rotation  stage  was  in  the  plane 
of  the  grating  so  that  the  laser  spot  samjded  the  same  area 
of  the  grating  throu^iout  a  scan.  Grating  depths  were 
measured  in  cross  section  with  a  scanning  electron  micro¬ 
scope  (SEM).  This  introduces  some  errors,  estimated  at 
±5  nm,  owing  to  uncertainties  in  the  SEM  calibration 
and  variations  in  the  grating  depth  for  the  different  areas 
sampled  in  the  optical  and  SEM  measurements. 
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Fif.  2.  Zero-order  reflectance  at  633  nm  for  610-nm-period 
gratinga  with  varying  grating  depth*.  The  left  side  [(a)-(d)]  pre¬ 
sents  experiment^  results;  the  right  side  ((e)-(h)]  presents  theo¬ 
retics]  modeling.  See  the  text  for  detaib. 


(3)  The  shift  in  the  SPW  resonance  to  smaller  angles 
with  increasing  depth, 

(4)  The  residua]  SPW  coupling  even  at  large  grating 
depth/grating  period  ratios, 

(5)  The  increasing  coupling  to  the  n  =  -1  diffraction 
order.  (The  sharp  spike  at  ~38*  corresponds  to  the  coUec- 
tion  of  the  n  s  -1  diffraction  order  in  the  optical  system 
and  is  not  part  of  the  zero-order  reflectance.  It  does  pro¬ 
vide  a  useful  monitor  of  the  energy  in  the  diffracted  order.) 

SEM  images  of  the  gratings  used  for  the  measurements 
of  Fig.  3,  taken  in  cross  section,  are  shown  in  Fig.  4. 
Note  that  the  grating  shapes  are  sinusoidal  at  low  depths 
but  gradually  show  increasing  harmonic  components  and 
trend  toward  rectangular  profUes  for  the  deepest  gratings 
investigated.  This  profile  modification  results  from  the 
grating  fabrication  technique.*^ 

Similar  experiments  were  performed  for  gratings  with 
periods  of  392  and  440  nm,  in  an  investigation  of  the  de¬ 
pendence  of  the  coupling  efficiency  on  grating  depth  and 
grating  period.  All  the  experiments  were  carried  out  at  a 
wavelength  of  633  nm,  in  avoidance  of  variations  in  the 
metal  optical  properties.  Qualitative^r  similar  behavior 
was  observed,  with  comparable  colliding  efficiencies  oc¬ 
curring  at  shallower  grating  depths  for  liner  gratings. 
Experimentally  determined  coupling  efficiencies,  reso¬ 
nance  angles,  and  resonance  linewidths  for  all  three  sets 
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EXPERIMENTAL  RESULTS 
Zero-order  reflectance  measurements  for  a  series  of  510- 
nm  gratings  with  increasing  depth  are  shown  on  the  left- 
hand  side  of  Fig.  2.  The  theoretical  modeling  shown  on 
the  right-hand  side  will  be  discussed  below.  The  mi^r 
features  to  be  noted  in  these  measurements  include  the 
following: 

(1)  The  excitation  of  SPWs  at  9  -  11.6*  cwresponding 
to  the  sharp  dip  in  the  reflectivity  (note  that  there  is  no 
diffracted  order  at  this  anide,  so  that  this  decrease  axm- 
qxHids  directly  to  energy  coiqited  into  the  SPW  modeX 

(2)  The  rapid  increase  in  coupling  efficiency  with  in¬ 
creasing  grating  depth  to  a  maximum  observed  coiqiling 
of  94%  at  a  grating  depth  of  35  nm, 

(3)  The  horizon  for  the  n  =  -1  diffracted  order  at 
$  —  13.8*  (this  is  apparent  as  the  cusp  in  the  reflectivity 
as  energy  is  transferred  from  the  specularly  reflected 
beam  to  the  diffracted  beam). 

Similar  results  for  deeper  gratings  are  presented  on  the 
left-hand  side  of  Fig.  3.  Note  the  lan^  angular  scale 
in  this  Hgure.  Specific  features  to  be  noted  include  the 
following: 


INCIDENT  ANGLE  (deg) 

Fig.  3.  Cmtinuation  of  Fig.  2  tar  deeper  grating*.  Note  the  ex- 


(1)  The  relatively  gradual  decrease  in  the  coiqding  effi¬ 
ciency  of  SPW’*, 

(2)  The  clear  broadening  of  the  SPW  resonance  with  in¬ 
creasing  depth. 


pended  angular  scale.  Again,  the  left  side  [(a)-(d)]  present*  ex¬ 
perimental  result*;  the  ri^t  *ide  [(e)-(h)L  theoretic^  modeling. 
Toe  sharp  spikes  on  the  experimental  data  at  —38*  correqxmd  to 
the  n  -1  difrractkm  ordw  entering  the  collection  optics  and 
are  not  a  part  of  the  zeru-order  reflectivity. 


Fig.  S.  Coupling  ofneiency  into  the  SFW  mode  for  gntingi  with 
pwioda  of  (a)  3n  nm,  (b)  440  nm,  and  (e)  510  nm  and  varying 
deptfaa.  Tte  eolid  curvea  are  the  reauh  of  thewetical  modeling; 
aee  the  text  for  detaDa. 
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Fig.  7.  Same  aa  Pig.  6;  except  here  for  SPW  reemanee  linewidth 
(FWHM). 
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GRATING  DEPTH  /  GRATING  PERIOD 
Fig.  8.  Data  of  Figs.  5-7  for  the  SPW  resonance  angle,  the  reso¬ 
nance  linewidth,  and  the  coupling  efficiency  refdotted  against  the 
dimensionless  parameter  h/d.  Note  that  this  parameter  is  ap> 
ivoximatefy  invariant  for  these  three  grating  periods. 


variations  of  grating  depth  and  {unfile  must  be  investi¬ 
gated  for  optimization  of  the  polarization  of  the  reflected 
beam.  Such  polarization  behavior  for  square  gratings 
was  predicted*^  and  demonstrated  in  the  infrared*’;  to 
our  knowledge,  though,  this  is  the  flrst  observation  of 
these  effects  in  the  visible  spectral  region.  SEM  pictures 
of  the  gratings  are  shown  in  Fig.  10. 

Further  increases  in  grating  depth  result  in  rectangu¬ 
lar  profiles  with  decreasing  line/space  ratios.  Angular  re¬ 
flectance  scans  for  deep  rectangular  gratings  are  shown 
in  Fig.  11  [(a)  320-nm  depth,  392-nm  period;  (b)  330-nm 
depth,  440-nm  period;  (c)  300-nm  depth,  510-nm  period]. 
There  is  a  broad  absorption  of  TM-polarized  ra^ation, 
while  a  large  diffraction  efficiency,  increasing  with  in¬ 
creasing  periods,  is  observed  for  TE  polarization.  Fig¬ 
ure  12  shows  SEM  pictures  of  these  gratings. 

THEORY 

During  the  past  thirty  years  many  grating  theories  based 
on  the  vector  character  of  the  electromagnetic  field  have 
been  developed.  An  excellent  summary  can  be  found  in 
Refs.  6, 13,  and  20.  Integral  methods,  developed  by  Petit 
and  others,  evaluate  the  field  at  any  point  in  terms  of 
an  integral  over  the  grating  surface  A  (cf.  Fig.  13).  A 
variation  of  these  integral  techniques  based  on  a  Green’s- 
function  formalism  specifically  directed  to  the  evaluation 
of  the  SPW  dispersion  relation  on  a  periodic  surface  was 
extensively  developed  by  Mills  and  co-workers^”*’  and 
Otagawa.”’  In  addition,  a  differential  formalism,  devel¬ 
oped  by  Nevierre  ei  al.,^  was  applied  to  SPWs  on  a  peri¬ 
odic  surface  by  Numata.””  All  Uiese  approaches  are  fully 
rigorous,  hold  for  arbitrary  profiles,  and  require  exten- 


of  gratings  are  presented  in  Figs.  5,  6,  and  7,  respectively. 
The  solid  curves  are  theoretical  and  will  be  discussed 
below.  With  increasing  grating  depth  each  ofthe  data  sets 
displays  an  initial  rapid  increase  in  coiqiling  efficiency, 
peaking  at  over  909h  and  a  slower  decrease  in  efficiency; 
an  approximately  quadratic  decrease  in  the  resonant  cou- 
fding  angle;  and  an  ai^roximately  quadratic  increase  in 
the  resonance  width. 

These  results  are  summarized  in  Fig.  8,  which  shows  all 
three  sets  of  data  {dotted  against  the  dimensionless  {»- 
rameter  h/d,  i.e.,  grating  depth/grating  {leriod.  Wit^ 
experimental  uncertainties  these  results  appear  to  follow 
a  common  bdiavior.  Some  of  the  variability  may  well 
arise  from  differing  grating  {nrofUes,  es|>ecially  fnr  the 
deeper  gratings. 

Further  increases  in  grating  depth,  accompanied  fay  a 
fhange  in  {anfile  to  rectangular,  le^  to  an  elimination  of 
SPW  effects.  For  approximately  square  gratings  strong 
pdarizatkm  effects  demonstrate  the  jiossibility  of  the  fab¬ 
rication  of  reflective  (wlarizers  for  visible  radiation.  This 
is  shown  in  Fig.  %  ahere  reflectance  scans  for  the  three 
grating  {leriods  are  given  f<v  depths  of  200  nm  (392-nm 
iwriod),  170  nm  (440-nm  period),  and  190  nm  (610-nm 
(leriod),  which  remihed  in  maximum  polarization  effects. 
Note  t^  for  TM  {xdarization,  almost  100%  of  the  inci¬ 
dent  miergy  (->96%  for  the  510-nm-period  grating)  is  cou¬ 
pled  into  the  first-order  diffractimi  i>eak  for  arq^  beyond 
the  horizon  for  this  order.  In  contrast,  only  ap{>roxi- 
mately  20%  of  the  energy  {xdarized  in  the  TE  dire^km  is 
coufded  out  of  the  sero-m^r  reflected  beam.  Careful 
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Fig.  ft  Sequence  of  zero-order  reflectance  scans  for  approxi¬ 
mately  square  gratings  (h/d  —  0.6)  showing  strong  polvizing 
effects:  (a)  A  *  200  nm,  d  «  392  nm;  (b)  k  s  170  nm,  d  > 
440  nm;  and  (c)  A  «  190  nm,  d  »  510  nm. 
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aive  computation,  typically  involving  matricea  of  order  40 
or  higher. 

A  number  of  differential  ana^yaea,  uaually  baaed  (m  a 
perturbation  expanaion  in  k^h,  where  h  ia  the  grating 
depth,  and  relying  on  a  Riqrleigh  or  plane-wave  annooch 
have  also  been  presented.****^  In  these  treatments  the 
scattered  fleld  anqditudes  are  alao  treated  as  smaQ  quan¬ 
tities,  of  the  order  of  k^h  times  the  incident  and  reflected/ 
transmitted  Held  amplitudes.  These  models  have  the 
virtue  of  relatively  sin^  ana^k  results  and  reaily  physi¬ 
cal  interpretation.  These  perturbational  approaches 
clearly  break  down  as  the  coupling  effkiencies  into  SPW 
modes  and  diffraction  orders  approach  unity,  since  they 
do  not  self-consistentfy  describe  the  necessary  decrease 
in  the  sero-order  reflected  and  transmitted  beams. 
Yamashita  and  Tsiyi**  treated  the  coupling  problem 
within  the  Rsyleigh  hypothesis  without  making  tite  smaD- 
signal  approximation  for  the  resonant  scattered  fields; 


however,  their  treatment  retained  the  expansion  in  k^h, 
adiich  limits  its  applicability  to  larger  grating  amplitudes. 
Nevertheless,  their  model  provided  an  elegant  analytic 
result  that  included  the  quadratk  [«(koA)*]  increase  in  the 
SPW  intensity  with  grating  dq>th  for  shallow  gratings, 
the  saturation  of  the  SPW  intensity  at  a  co\q>ling  efH- 
ciency  near  100%,  and  a  gradual  decrease  of  the  coiq>ling 
{«(koh)~*]  for  deeper  gratings  as  the  energy  is  coupled 
back  into  the  radiating  fields,  in  qualitative  agreement 
with  the  experimental  results  presented  above,  l^ber^ 
presented  a  related  analysis,  based  on  the  extinction  theo¬ 
rem,  that  retains  the  nonresonant  terms  in  a  perturbation- 
theory  expansion  in  koh  while  also  keeping  the  resonant 
terms. 

These  models  begin  with  a  time-harmonic,  plane-wave 
expansion  of  the  electromagnetic  field  in  the  regions  out¬ 
side  the  grating  kerf  (z  <  0  and  z  >  h.  Fig.  13);  the  mag¬ 
netic  Helds  of  the  TM-polarized  Helds  are  given  fay 

B*  =  B'e,  =  {exp[i{*,y  +  *.z)] 

+  2  B",  exp(i(*,y  -  a,z)]}e„  z  <  0  (3a) 


and 


B"  =  B"ex  =  2  exp[i(*.y  +  /5,z)]e„ 


z  >  h,  (3b) 

where  k,  -  kp  sin  B,  k,  =  kp  cos  B,  k,  =  k,  +  ng, 

n  =  0,±1,  ±2 . g  =  2v/d  is  the  grating  wave  vector 

with  d  the  grating  period,  a.  =  i(k,*  -  kp*)‘'*,  = 


Incident  Angle  (deg.) 

Fig-  11-  Sequence  ot  zero-order  reflectance  scanf  for  deeper 
rectangular  gratings  {h/d  -  t)  showing  strmg  absorption  effects: 
(a)  h  *  Si20  nm,  d  »  3K  nm;  (b)  A  »  330  nm,  d  «  440  nm;  and 
(c)  A  >  300  nm,  d  ■>  610  nm. 
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Fig.  12.  SEM  pictUTM  of  tlw  gntingi  OMd  for  the  nMuuie- 
ments  of  Fig.  IL 


and  Cadilhac**  showed  that  for  TE  polarization  the  Rqr- 
leigh  expansion  is  convergent  for  h/d  <  0.14,  and  pointed 
out  that  reliable  results  could  be  obtained  even  tor  h/d  val¬ 
ues  as  great  as  twice  this  limit.  A  detailed  analysis  of  the 
applicability  of  the  Rayleigh  hypothesis  to  real  materials 
was  recently  presented.”  Glass  et  al.**  pointed  out  that 
the  extinction  theorem  formalism  leads  to  the  identical 
dispersion  relaticm. 

'the  i4>proach  used  here  in  the  comparison  with  experi¬ 
ment  is  to  apply  the  Rayleigh  hypothesis  for  TM-polarized 
input  radiation.  For  the  gratings  used  in  these  experi¬ 
ments,  h/d  extends  to  0.8,  althou^  detailed  coD[^>ari8on8 
are  attempted  only  to  —0.3.  The  calculation  does  not 
make  a  small-signal  approximation  for  either  the  grating 
depth  or  the  diffracted  field  amplitudes.  The  plane-wave 
expansion  [Eqs.  (3)]  is  truncated  fay  the  retention  of  reso¬ 
nant  terms,  e.g.,  n  =  0  and  n  =  -1.  In  addition,  the  next 
terms,  n  =  -2  and  n  =  +1  are  retained  in  the  numerical 
evaluations  and  are  determined  to  be  small  relative  to  the 
resonant  terms.  Energy  conservation,  i.e.,  a  constant 
total  of  the  energies  in  the  diffracted  beams  and  absorbed 
in  the  metal,  is  used  as  a  further  check  on  the  calculation. 

Thus  Eq.  (3)  is  assumed  to  hold  up  to  the  grating  sur¬ 
face  defined  fay 

f(y)  =  u  sin(gy),  (4) 

i^ereu  =  A/2  is  the  grating  anqilitude.  Using  the  gener¬ 
ating  function  for  Bessel  functions,  we  can  write 

e^viiyfiy)]  =  sin(gy)]  «  2  exp(ipgy)J,(yu) . 

(5) 

The  boundary  conditions  satisfied  these  fields  are 

=  (6) 


dB’iy.z)  1  dB-(yiz) 
bn  c.  bn 

where 


Fig.  13.  Geometry  used  in  the  analysis. 


i(A,*  -  sml  Cm  is  the  metal  dielectric  constant. 

Within  the  grating  kerf  (0  <  z  <  A),  the  validi^  of  this 
expansion  is  not  weU  established.  Rayleigh  made  the 
assunqitiim,  known  as  the  Rasdei^  hypothesis,  that  for 
sufficiently  shallow  gratings  this  expansion  is  valid  every¬ 
where.  This  hypothesis  sraa  investigated  fay  Vhn  den  Berg 
and  Fokkema”  and  shown  to  be  analytically  correct;  Petit 


Applying  these  boundary  conditums  leads  to  an  infinite 
set  of  couided  linear  equidions: 

Jp(k,u)  +  =  0 

(8) 

and 

(a.  -  -  2[«.  -  V,-. 

X  («.«)  -  i  -  0, 

(9) 


sHiere  p  and  n  are  integers  extending  fhnn  to  +«. 
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Fig.  14.  Comparison  of  (a)  cakulationa  of  Garcia’*  with  ((b)-(e)] 
present  calculations  for  800-nm-period  Au  gratings  at  633  nm 
(A,  =  h/2d). 


Note  that  no  small-gignal  a^qjroximation  was  made  in  the 
derivation  of  these  equations.  These  equations  are  now 
truncated  by  the  retention  of  only  the  fields  for  ft  =  0  and 
n  =  -1,  the  reflected/transmitted  and  resonant  diffracted 
terms,  which  are  expected  to  be  large  based  on  [diase- 
matz-hing  arguments  [cf  Eq.  (2)  and  Fig.  1].  In  addition, 
the  fields  for  n  =  4-1  and  n  s  -2  were  retained  as  a  check 
on  the  convergence.  The  relative  Held  intensities  for 
these  two  orders  for  a  SlO-nm-period  grating  at  a  d^th  of 
22  nm  were  0.017  and  0.0018,  respectively,  lending  support 
to  this  truncation  procedure.  This  leads  to  solving  an 
8x8  matrix  inversion  for  the  field  intensities. 

Using  the  Green’s-fiinction  approach  developed  by  Toigo 
c(  of.,**  Garda  evaluated  the  Sl^  fields  and  1^  shsqw  fw 
frimifniHal  gratings.**  The  results,  obtained  by  the  inchi* 
sion  of  60  terms  in  the  numerical  analysis,  are  shown  in 
Fig.  14(a)  foe  parameters  apinopriate  to  an  800-nm-period 
Au  grating  of  various  depths  h,  =  A/2d.  For  this  grating 
period  and  wavelength  (633  nm),  X/i  <  1,  and  there  are 
two  possible  propagating  diffraction  orders.  The  reso¬ 
nance  shown  in  Fig.  14  corresponds  to  the  n  =  -f  1  SPW 
coupling;  there  is  also  a  propagating  n  -1  diffraction 
order  throu^iout  this  angular  range.  The  numbers  label¬ 
ing  the  SPW  curves  in  the  figure  represent  the  relative 
intensity  of  the  mode.  The  results  frmn  the  present  treat- 
mmit  fm:  the  same  parameters  are  shown  in  Figs.  14(b)- 
14(d).  Note  the  overall  similarity  of  the  calculated 
intensities  and  line  shapes.  Interestingly,  the  largest 
discrepancy  is  for  the  shallowest  grating,  ht  -  0.01,  wfame 


the  present  calculation  yields  a  SPW  intensity  that  is 
24%  larger  than  the  more  rigorous  calculation  of  Garcia.** 
For  deeper  gratings  the  agreement  is  generally  within 
10%.  A  similar  result  was  found  by  Weber**  in  his  devel¬ 
opment  of  a  coupled-mode  anal^is.  This  comparison 
indicates  that  despite  the  relative  simplicity  of  the  formu¬ 
lation  the  present  model  provides  sul»tantial  insight  into 
SPW  coupling. 

COMPARISON  OF  THEORY 
AND  EXPERIMENT 

The  initial  step  in  a  comparison  of  this  model  with  the  ex¬ 
perimental  results  is  the  establishment  of  the  deposition- 
dependent  dielectric  properties  of  the  Ag  films.  Various 
values  of  the  Ag  dielectric  constant  at  a  633-nm  wave¬ 
length  have  been  reported,  viz.,  (-16.4,0.54),’* 
(-17.9^0.49),*®  and  (-18.91 2.35)."  As  is  reflect^  in  this 
wide  dispersion  of  reported  values,  the  dielectric  prop¬ 
erties  of  thin  films  are  dependent  on  deposition-  and 
substrate-dependent  columnar  structure,  granularity, 
subideal  density,  and  incorporated  inqnirities.  Given  this 
uncertainty,  the  approach  adopted  here  is  to  acljust  the 
dielectric  properties  to  obtain  a  best  fit  to  the  experimen¬ 
tal  result  for  a  single  grating  [510-nm  period,  22-nm  depth; 
Fig.  2(c)],  which  gives  a  value  of  (-13.65,0.82);  this  value 
was  used  for  all  the  comparisons  between  theory  and 
experiment.  The  fit  for  this  grating  is  shown  in  Fig.  15; 
excellent  agreement  is  obtained  for  a  calculated  depth  of 
19  nm.  In  subsequent  measurements,  which  will  be  re¬ 
ported  in  detail  elsewhere,  the  dielectric  properties  of  Ag 
films  deposited  on  Si  substrates  under  identical  conditions 
were  determined  independently  by  ellipsometric  tech¬ 
niques  to  be  (-12.9(0.76),  within  7%  of  the  values  obtained 
here.  The  grating  fabrication  and  the  film  deposition  in¬ 
volved  in  the  present  8tuii(y  were  carried  out  over  a  period 
of  several  months,  and  deposition  conditions  and  ^  re¬ 
sulting  film  properties  invariably  drift  over  this  time 
period;  no  abetment  was  made  in  the  dielectric  values 
used  in  the  comparison  between  theory  and  experiment. 


INCIDENT  ANGLE  (deg) 

Fig.  1&  Calculated  zero-order  reflectivity  SPW  resonance  line 
shapes  for  c  -  (-1386,0.82)  at  grating  depth  A  =  19  nm.  The 
expoimental  rnult  (h  =  22  nm)  is  shown  as  the  solid  curve. 
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Fig.  16.  Calculated  SPW  reaonance  parameters  platted  against 
the  dimensionless  parameter  h/d. 


Comparisons  with  the  experimental  resiilts  are  shown 
on  the  right-hand  sides  of  Figs.  2  and  3.  In  each  case  the 
grating  depth  for  the  calculation  is  adjusted  to  give  cou¬ 
pling  efficiencies  that  match  the  experimental  results. 
There  is  overall  excellent  agreement.  For  the  deepest 
grating  shown  experimentally,  A  =  145  nm,  there  is  a 
substantial  difference  with  the  model  that  is  probably 
associated  with  the  significant  deviation  from  a  simple 
sinusoidal  profile  at  this  depth  (cf.  Fig.  4).  Figures 
present  the  model  results  along  with  the  experimental 
data  for  the  coupling  efficiency,  the  resonance  an^de,  and 
the  resonance  linewidth  fw  aQ  three  grating  periods  invM- 
tigated.  There  is  an  excellent  qualitative  match  between 
theory  and  experiment  The  most  significant  deviatka, 
that  for  the  linewidth  of  the  440-nm-period  gratings,  is 
most  probably  due  to  a  variation  in  the  film  dielectric 
properties  fw  this  set  of  gratinga  The  couiding  efficien¬ 
cies  first  increase  rapidly  as  the  grating  depth  is  increased, 
saturate  near  100%  coupling,  and  then  decrease  as  the 
resonance  becomes  overcoupled  owing  to  radiative  damp¬ 
ing  back  into  the  zero-order  radiation  fields.  This  cou- 
j^ng  change  is  accompanied  fay  aiq>roximately  qu^ratic 
changes  in  the  resonance  an^e  and  linewidth  with  in¬ 
creasing  grating  depth.  Remarkably  similar  qualitative 
behavior  was  predicted  Iqr  Weber  and  Mills**  in  their 
numerical  ana^is  of  grating  coupling  for  800-nm-period 
Ag  sawtooth  gratings.  Thqr  evaluated  the  third-order 


coupling  at  A  =  480  nm  and  found  a  maximum  coupling 
strength  for  grating  depths  of  40-60  nm. 

The  model  results  are  summarized  in  Fig.  16,  vdiich 
shows  the  calculated  resonance  angle,  resonance  line- 
width,  and  cotq>ling  strei^lth  plotted  against  the  dimen¬ 
sionless  parameter  h/d  for  the  three  grating  periods 
investigated.  For  evaluation  of  the  coupling  efficiency 
this  parameter  is  reasonably  invariant;  there  are  more 
significant  deviations  in  the  evaluation  of  the  resonance 
line-shape  parameters. 

Overall,  this  simple  model  provides  a  good  picture  of  the 
experimentally  observed  resonance  variations.  Disagree¬ 
ments  between  theory  and  experiment  increase  for  deeper 
gratings  and  larger  periods.  A  significant  phenomenon 
not  included  in  the  theoretical  model  is  the  deviation  from 
sinusoidal  grating  profiles,  which  increases  as  the  grating 
depth  is  increased  and  also  is  more  significant  for  larger 
grating  periods.  The  Hlms  also  exhibit  significant  sur¬ 
face  roughness  (see  Fig.  4),  which  has  not  been  included  in 
the  model  and  may  influence  the  observed  line  shapes. 

Calculated  zero-order,  TM-polarization  reflectivity 
curves  for  deeper  gratings,  up  to  h/d  ~  1,  are  shown 
in  Fig.  17.  Again,  there  is  good  qualitative  agreement 
with  the  large  coupling  into  the  diffraction  order  for 
deep  gratings  (h/d  ~  0.5),  although  detailed  comparisons 
are  not  possible  because  of  the  strong  deviation  from  a 
simple  sinusoidal  profile  of  the  experimental  gratings. 
The  model  does  not  show  the  absorption  of  TM-polarized 
radiation  seen  experimentally.  A  detailed  comparison 
with  experiment  requires  deep  sinusoidal  gratings,  which 
are  inherently  difTicult  to  fabricate  because  of  the  strong 
nonlinearities  of  existing  photoresists.**  More  nearly 
sinusoidal  gratings  can  be  fabricated  on  a  transparent 
substrate;  research  continues  in  the  pursuit  of  a  better 
experimental  test  of  the  model  fw  deep  gratings.  The 
theoretical  model  loses  much  of  its  simplicity  for  rectan¬ 
gular  gratings,  where  many  Fourier  conqwnents  of  the 
grating  profile  are  comparable  in  intensity. 


INCIDENT  ANGLE  (deg) 

Pig.  17.  Calculated  zero-order  reflectivity  for  deep  (h/d  ~  1)  ri- 
nuaoidal  gratings  in  TM  pcdarisation  (d  =  610  nm). 
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SUMMARY 

A  systematic  experimental  and  theoretical  8tu<i|y  of  first- 
wder  grating  coupling  to  SPW’s  existing  at  an  air-Ag 
interface  was  reported.  The  experiment  extends  the 
range  of  investigated  grating  depths  to  h/d  ~  1.  For  TM- 
polarized  incident  radiation  this  includes  the  entire  range 
of  SPW-radiation  coupling— from  underdamped,  to  nearly 
100%  coupli>'ig,  to  overdamped  and  the  ultimate  di8aH>ear- 
ance  of  the  SPW  resonance  from  the  observed  zero-order 
reflectance.  Strong  polarization  and  absorption  effects 
are  observed  for  the  deepest  gratings.  A  simple  theoreti¬ 
cal  model,  based  on  the  Rayleigh  hypothesis  and  retaining 
cmly  resonant  diffraction  terms  without  making  a  small- 
signal  ai^roximation,  provides  good  agreement  with  the 
experimental  results. 

Several  extensions  of  this  work  are  immediately  af^r- 
ent.  For  larger  grating  periods  there  are  more  di^rac- 
tion  orders  and  SPW  coupling  resonances.  Interesting 
coupling  effects  occiu:  when  two  of  these  resonances  occur 
at  ai^roximately  the  same  angle.  These  effects  have 
been  investigated  experimentally  and  theoreticaUy  and 
win  be  reported  in  a  subsequent  publication.  In  addition, 
the  SPW  resonance  can  be  used  to  characterize  metal  op¬ 
tical  constants  under  a  variety  of  deposition  conditions. 
A  study  is  under  way  that  conqMues  this  technique  with 
more  conventional  techniques  such  as  eUipsometry. 
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For  grating  periods  larger  than  the  excitation  wavelength,  multiide-grating  orders  couple  incident  optical  radia¬ 
tion  to  the  surface  jdas^  waves  (SPW*s)  characteristic  of  the  metal-air  interface.  Fm  a  grating  period  that  is 
an  inters]  multiple  of  the  wave  vector  of  these  surface  modes,  two  resonatrces  become  d^nerate  in  co«q>ling 
angle.  There  are  also  permitted  diffraction  orders  at  this  coupling  an^.  The  vicinity  of  this  multiple-mode 
coupling  resonance,  sdrm  several  free-qiace  electromagnetic  modes,  as  well  as  two  surface  modes,  are  cotqded 
by  different  wders  of  a  grating,  is  known  as  a  minigap  region.  Not  surprisingly,  the  response  surface  displays 
complex  dependences  on  frequency,  angle,  and  grating  proHle.  A  detailed  experimental  and  thewetica]  study  is 
presented  of  the  optical  response  at  633  nm  in  the  (-t-l,  -2)  minigap  region  for  Ag  films  deposited  on  phoUditho- 
graphkally  defin^  870-nm-period  gratings.  Measurements  of  the  0-<nder  reflectance  and  the  -1-order 
Effraction  are  presented  for  a  wide  progression  of  grating  depths.  The  SPW  rescmances  depend  on  the  grating 
depth,  and  this  variation  is  used  to  tune  through  the  minigap  region  for  a  fixed  aravelength  and  period.  Simi¬ 
lar  measurements  are  presented  for  a  sin^  grating  as  a  function  of  wavelength  through  the  minigap  region. 
In  both  measurements  the  0-order  response  shom  only  a  single  broad  minimum  as  the  resonances  approach 
degeneracy,  arfaile  the  -1-order  diffractikm  shows  clearly  defined  momentum  gaps.  A  simple  theoretic^  model 
ba^  on  tite  Rayleiidt  hypothesis  is  presented  that  gives  a  good  qualitative  picture  of  the  response.  The  re¬ 
sponse  surfaces  are  sensitive  to  the  grating  proflK  and  detailed  modeling  requires  inclusion  of  higher-or^ 
grating  components. 


1.  INTRODUCTION 

Optical  interactions  with  metallic  gratings  have  long  fas¬ 
cinated  scientists.  Wood’s  initial  observations  of  grating 
anomalies'  were  classified  by  Fano*  into  (1)  diffraction 
anomalies  associated  with  the  redistribution  of  energy 
among  diffraction  orders  when  one  order  passes  over  the 
horizon  (i.e.,  the  propagation  direction  readies  an  an^e  of 
90*  to  the  grating  normal)  and  (2)  resonance  anomalies 
associated  with  the  propagation  of  guided  waves  akmg  the 
grating  interface.  A  large  botly  of  experimental  and  theo¬ 
retical  research  relating  to  optical  interactions  with  peri¬ 
odic  surfaces  has  been  accumulated.*-*  This  paper  is 
cmicemed  with  resonarKs  anomalies  and  is  focused  on  the 
propagation  of  guided  waves,  known  as  surface  {dasma 
wavea  (SPW’s),  on  bare  metallic  surfaces. 

For  grating  periods  leas  than  the  ezdtation  wavdmigth, 
on)y  fint-«tler  ooi^ng  to  SPW’s  is  pen.nitted  and  has 
bsMi  described  in  detail  elsewhere.*  For  grating  periods 
larger  than  the  incident  wavelength  it  is  possible  to  couple 
to  SPW’s  in  higfa«r  orders  and  to  obeo^  electromagn^ic 
interactions  between  these  coiqilings.  These  interactions 
are  the  ru^  focus  ct  this  pqier.  The  first  experimental 
evidmee  of  these  intmactioos  between  diflermt  coupling 
orders  was  provided  by  Stewart  and  Gallaway,*  who  ob¬ 
served  that  in  some  cases  Wood’s  anomalies  did  not  become 
coincident  but  replied  each  other  as  the  optical  wave- 
ler^th  was  varied  to  bring  the  resonances  closer  together. 
Ritdiie  c(  a/.*  also  obewed  g^M  in  the  dispoaion  relatkm 
of  the  SPW’s  and  attyibuted  them  to  hi^ier-order  interac¬ 
tions  between  SPW’s.  Hutley*  and  Hutley  and  Bird*  car- 
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ried  out  a  detailed  experimental  stu4y  of  the  anomalies  of 
sinusoidal  diffraction  gratings  in  the  visible  spectrum 
The  anomalies  were  characterized  in  a  Littrow  configura¬ 
tion  as  a  functicm  of  groove  depth,  period,  and  wavelength, 
and  in  some  cases  a  reluctance  of  anomalies  to  merge  was 
observed.  Kroo  et  ai.’^  observed  similar  gaps  in  the  dis¬ 
persion  relation  of  SPW’s  in  metal-oxide-metal  structures 
and  attributed  them  to  the  crossing  of  (+1,  -1)  coiqrling  to 
SPW’s.  Chen  et  a/."  experimentally  investigated  grating 
cou|ding  to  SPW’s  in  the  minigap  regions,  i.e.,  degenerate 
coigiling  to  SPW’s  in  two  different  grating  or^rs,  such  as 
(■fl,  -2)  or  (-t-2,  -2),  and  measured  energy  gaps  for  vari¬ 
ous  minigap  regiais.  Wsber  and  Mills'*  and  Tran  et  al.'* 
showed  throretically  that  interactkms  of  elastically  scat¬ 
tered  SPW’s  result  in  conqilez  reqxmse  surfaces  in  the 
vicinity  of  the  minigaps  and  that  a  direct  interpretation  in 
terms  of  a  modal  diipeieion  rdathm  is  difficult  at  best. 
Detailed  modetiiy  and  carefiil  evaluation  of  coufdings  in 
both  amplitude  and  phase  are  necessary  for  a  full  aipred- 
ation  of  the  experimental  results.  Heitman  et  al.'* 
showed  from  experimental  studies  on  metal-oxide-metal 
structures  that  the  existence  of  energy  or  momentum 
gaps  was  dependent  on  the  iiKKie  of  excitation  and  the  ex¬ 
perimental  observable:  i.e.,  using  fast  electton  excitation 
and  light  emitted  from  metal-oxide-metal  structures, 
they  observed  energy  and  momentum  gaps,  but  only  en¬ 
ergy  gaps  were  obsm-ved  from  optical  excitation  (similar 
to  Chen’s  results).  They  also  pointed  out  that  the  anoma- 
Ues  strongly  depend  on  grating  period,  dqiUi,  and  profile 
as  wril  as  on  the  excitaticm  mechanism  and  the  method 
of  detection. 
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Pic-  DUpenioa  reUtion  of  SPW’t  for  a  IomIom,  free-etectroo 
metal  (c  >■  1  -  ^  ***  normalued  to  at,  aad  k,  « 

mpfe.  Thefratiag  vectore  oorreapandiogto  the  a  «  ±1  and  A  * 
±3  ordera  of  a  grating  of  period  d  are  ihown  as  vertical  dashed 
lines.  The  range  of  wave  vectors  accessible  hy  varying  the  inp^ 
an^  firomnewmal  to  grating  incidence  for  n  =  -t-landn  »  -3ia 
shown  as  horizontal  lines.  Note  that  for  this  choice  of  parame* 
tars  (A/tf  -  0.73),  there  are  two  SPW  resonances  and  that  the  -1 
diffraction  or^  propagatea  throughout  the  angular  range. 

When  higher-order  couplings  to  SPW’s  are  studied,  the 
permitted  diffraction  orders  are  strongly  coupled  to  the 
SPW’s;  therefore  the  incident  beam,  all  propagating  dif¬ 
fracted  beams,  and  the  SEWs  must  be  ta^n  into  account 
for  a  complete  understanding  of  these  complex  interac¬ 
tions.  The  effect  of  SPW  resonance  on  diffraction  orders 
has  been  studied  by  many  authors.**~*^  However,  all  those 
etudies  dealt  with  sin|^  SPW  resonance;  the  effects  of 
multiple  couplings  have  not  been  investigated  experimen¬ 
tally.  Here,  a  aystematk  theoretical  and  experimental  in¬ 
vestigation  fd  SPW  interference  in  the  (4-1,  -2)  minigap 
region  is  presmited.  Both  the  0-order  reflectance  and  the 
-l-ordor  diffraction  (the  only  permitted  diffraction  order 
in  the  angular  region  of  SPW  coupling)  were  measured  at 
6S3  nm  aa  a  ftinetion  of  the  inddeot  angle  for  varioas  grat¬ 
ing  depths  and  profUee  at  a  fixed  grating  period  of  870  nin. 
The  angular  shifts  in  the  SPW  resonance  angles  as  a 
result  of  the  increasing  grating  depths  were  used  to  tune 
throu^  the  coupling  degeneracy.  By  adding  intenaitiea 
in  the  0-order  reflectance  and  the  -1-order  diffraction  to 
ftnd  the  total  reflected  energy,  we  obeerved  momentum 
gape.  Also,  for  a  sin^e-grating  profile  similar  nMasure- 
nMnts  were  repeated  as  a  ftinctkm  of  wavelength.  Con- 
aisteni  with  the  experimental  results  discussed  abovev 
momentum  gaps  were  not  found  in  the  0-ord«r  reflec¬ 
tance;  however,  wdl-defined  momentum  gape  were  ob¬ 
served  in  the  -1  diffraction  order. 

finally,  a  sisqdified  theoretical  formalism,  developed 
for  the  analyBis  of  fint^ader  coupling,' is  extwoded.  This 


analysis,  based  on  a  simple  truncation  of  the  plane-wave 
expansion  of  the  scattered  fields  within  the  Rayleigh 
hypothesis,  provided  remarkably  good  agreement  with 
the  first-order  coupling  results.  Not  surprisingly,  while  ^  ) 
the  mi^r  features  of  the  experiment  are  reproduced,  the 
agreement  is  lees  satisfactory  for  this  more  demand¬ 
ing  case. 

The  dispersion  relation  for  SPWa  at  a  metal-air  inter¬ 
face  is  given  by" 

*apw  =  +  <)]"*,  (1) 

where  ho  is  the  free-space  propagation  vector  and  t  is 
the  metal  dielectric  constant.  The  excitation  of  SPW’s 
takes  place  whenever  the  wave  vector  of  the  incident  light 
parallel  to  the  metal-air  grating  interface  satisfies  the 
condition 


ho  sin  0  -I-  2irn/d  -  ±kgem,  (2) 


vdiere  9  is  the  incident  angle,  d  is  the  grating  period,  and 
n  =  0,  ±1, ±2,...  represents  the  coupling  order.  For  a 
lossless  free-electron  metal  (e  =  1  -  the  disper¬ 

sion  relation  [Eq.  (1)]  is  plotted  in  Fig.  1,  where  k,  =  Wp/c. 
The  vertical  dashed  lines  represent  the  various  grating 
orders,  and  the  horizontal  lines  indicate  the  accessible 
wave  vectors  (2irn/(f  to  ho  4-  2vn/d)  in  the  4-1  and 
-2  ord -~8.  Note  that  coup'ing  to  SPW’s  can  be  achieved 
for  both  of  these  orders  for  the  choice  of  w  and  d  corre¬ 
sponding  to  the  present  experiment.  C(Hq>ling  in  the  1 
(-2)  orders  results  in  the  jHqpagatkm  of  SPW’s  in  the  for¬ 
ward  (backward)  [i.e.,  same  (opposite)]  direction  as  the  y 
component  of  the  incident  wave  vector.  The  -1  diffrac¬ 
tion  order  is  permitted  throu^iout  the  angular  range  of 
the  two  SPW  couplings,  and  its  properties  are  stnm^  in¬ 
fluenced  by  the  SPW  resonances.  Also  rote  that  degener¬ 
acy  in  the  coupling  angle  9  correqwnds  to  the  coupling  of 
the  oppositely  directed  SPW’s  by  an  integral  multiple  of 
the  grating  wave  vector  (3  in  this  ease).  In  a  fashion  fa¬ 
miliar  from  many  examples  of  mode  coupling,  this  degen¬ 
eracy  results  in  the  interactions  between  the  SPW  modes 
that  give  rise  to  the  minigap  effects. 


Fig.  3.  ExperimsBtalamngsimntfiarrcfketkaaiiddifftaetiao- 
ardsr  aasasursments. 
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Fi(.  3.  0-order  renectancee  at  633  run  for  870-iun-period  grat- 
inga  with  varying  grating  deptha  h.  The  left-hand  column 
preaenta  experiment^  reaulta;  the  i^t-hand  cohunn  preaenta 
theoretical  modeling.  See  the  text  for  detaila. 


2.  EXPERIMENT 

Gratings  were  fabricated  hologra{diical]y  in  positive  photo¬ 
resist  on  Si  substrates  with  a  single-frequency  Ar-ion 
laser  smarce  at  488  nm.**  Grating  depths  were  varied  by 
adjustment  of  the  exposure  time.  This  results  in  varia¬ 
tions  in  both  the  depth  and  the  grating  profile,  which 
varies  from  approximately  sinusoidal  at  small  depths  to 
rectangular  at  the  deepest  gratings  investigated.  These 
photoresist  gratings  were  coated  with  -lOO-nm-thick  Ag 
films  with  the  use  of  an  electron-beam  ev^aorator  system 
at  room  temperature  and  at  a  background  vacuum  rfbw- 
lO-*TDrr. 

The  optical  measurements  were  scans  of  reflectance 
(diffraction)  versus  incident  angle  in  the  experimental 
setup  shown  in  Pig.  2.  In  this  arrangement  two  concen¬ 
tric  rotation  stages  are  used,  with  the  sample  mounted  on 
the  inner  stage  and  the  detectw  mounted  on  the  outer 
stage.  The  data  acquisition  step  invcdved  a  personal  com¬ 
puter  interfaced  wi^  a  stepper  motor  controller,  which,  in 
turn,  cmitrolled  the  two  rotation  stages.  Most  measure¬ 
ments  were  made  at  the  He-^e  laser  wavelength.  For  the 

Ovariabte-wavdength  measurements  a  cw-dye-laser  source 
was  used.  Grating  dqiths  were  measured  in  cross  section 
with  a  scanning  electro  microscope  (SEM).  This  intro¬ 
duces  errors,  estimated  at  ±5  nm,  owing  to  uncertainUes 
in  the  SEM  calibration  and  to  variations  in  grating  depths 


for  the  different  areas  sampled  in  the  optical  and  SEM 
measurements. 

3.  O-ORDER  REFLECTANCE 
MEASUREMENTS 

Figures  3a-3d  show  the  O-order  reflectance  scans  as  the 
grating  depth  is  increased.  The  plots  in  Pigs.  3e-3h  are 
obtained  from  the  modeling  and  will  be  discussed  in  Sec¬ 
tion  8.  The  following  important  features  should  be  noted: 

(1)  The  excitation  of  SPW’s  in  the  first  (n  =  +1)  and 
second  (n  »=  —2)  orders  at  —18"  and  24.7",  respectively, 

(2)  The  horizon  for  the  associated  diffraction  orders  at 
16.8"  and  27.1",  respectively, 

(3)  The  increase  in  the  first-order  coupling  efficiency  as 
the  grating  depth  is  increased, 

(4)  The  relative  weakness  of  the  second-order  coupling, 
(6)  The  decrease  in  the  angular  separation  of  the  two 

coupling  resonances  as  the  grating  depth  is  increased. 

Also  note  that  the  angular  shifts  in  the  second-order  reso¬ 
nance  are  almost  twice  those  in  the  first  order. 

For  deeper  gratings  the  corresponding  measurements 
are  shown  in  Figs.  4a-4d;  as  before.  Figs.  4e-4h  are  ob¬ 
tained  from  theory.  Note  the  foDowing  features: 

(1)  The  increase  in  the  -2-order  coupling  strength, 

(2)  The  merging  of  the  first-  and  second-order  cou¬ 
plings,  resulting  in  a  saddle-point  line  shape. 


INCIDENT  ANGLE  (deg) 

Fig.  4.  Same  as  Fig.  3  but  for  deeper  gratings. 
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Fig.  &  CntM-Mction  SEM  photognpiu  of  the  gratingi  used 
for  the  esperimente.  The  meeaured  depths  ere  (a)  95  nm, 
(b)  110  nm,  (c)  122  nm,  and  (d)  150  nm. 


(3)  The  emergence  of  a  broad  minimum  (Fig.  4d), 
where  Tint-  and  second-order  couplings  can  no  longer  be 
distingiiished. 

This  is  similar  to  the  behavior  obsm^red**  and  calculated** 
for  similar  fured-frequeD<7  reflectivity  scans  in  the  vicin¬ 
ity  of  a  minigap.  As  a  re^t  of  the  saddle-point  response 
surface,  only  a  sin|d^  broad  resonance  is  observed  v^en  0 
is  varied  for  fixed  tt,  while  two  distinct  resonances  and  an 
energy  g^p  are  found  when  w  is  varied  for  fixed  A 
In  Fig.  6  SEM  (nofiles  of  the  gratings  used  to  obtain 
Figs.  4a-4d  are  shown.  Increasing  grating  depth  also 
results  in  increasing  deviation  from  sinusoidal  profiles  as 
a  result  of  the  grating  fabricatkm  technique.'*  For  stiQ 
deeper  gratings  the  corresponding  measurements  are 
shown  in  Fig.  6,  where  the  vertical  scale  is  unchanged 


from  that  of  Figs.  3  and  4.  For  these  approximately 
rectangular  grating  structures  (cf.  Fig.  7),  there  is  an 
overall  decrease  in  the  0-order  reflectance,  an  absence  of 
clear  SPW  coupling  resonances,  and  line  shapes  that  do 
not  lend  thems^es  to  simple  interpretation. 

4.  -1-DIFFRACTION-ORDER 
MEASUREMENTS 

For  this  grating  period  and  measurement  wavelength  the 
-1  diffraction  order  is  a  propagating  mode  throu^unit  the 
entire  incident  angle  range  including  that  of  the  SPW  reso¬ 
nances.  Coupling  to  this  propagating  order  provides  addi¬ 
tional  information  on  the  SPW  resonance  characteristics. 
In  fact,  the  first  observations  of  Wood’s  anomalies  were 
related  to  diffraction  orders  rather  than  to  0-order  reflec¬ 
tance.  Measurements  of  the  energy  coupled  into  the  SPW 
mode  must  take  into  account  both  the  reflectance  and  the 
energy  coupled  into  the  diffraction  order.  DiPTraction- 
order  scans  corresponding  to  the  0-order  scans  of  Figs.  3a- 
3d  are  shown  in  Figs.  8a-8d;  as  above,  the  plots  in 
Figs.  8e-8h  are  obtained  from  theory.  The  vertical  scale 
on  these  measurements  is  the  same  as  that  for  the  0  order, 
so  the  diffracted  energy  in  the  -1  order  is  significant 
when  compared  with  the  reflectivity.  The  following  fea¬ 
tures  should  be  noted; 

(1)  The  line  shapes  are  derivativelike  as  the  incident 
angle  is  varied  through  the  SPW  resonances, 

(2)  The  coupling  strengths  are  equal,  in  contrast  with  the 
much  weaker  -2-order  coupling  (Akserved  in  the  reflectiv¬ 
ity  (cf.  Figs.  3a-3d), 

(3)  Neither  the  maxima  nor  the  minima  in  the  diffraction- 
order  intensities  correspond  to  the  0-order  minima, 

(4)  There  is  an  enhancement  of  the  diffraction-order  in¬ 
tensity  in  the  region  between  the  SPW  interactions, 

(5)  As  the  grating  depths  increase,  the  angular  separa¬ 
tion  between  the  SPW  resonances  as  measured  between 


10  15  20  25  30 

INCIDENT  ANGLE  (deg) 

Fig.  6.  0-order  reflectances  at  633  nm  for  gratings  with  (a)  resid¬ 
ual  SPW  effects  and  (b)  absorptive  bdiavior. 
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(3)  The  enhancement  in  the  -1  diffraction  order  is  such 
that  the  minimum  in  the  0  order  (Fig.  4d)  corresponda  to 
the  maximum  in  the  diffraction  order  (Fig.  9d). 

A  comparison  of  the  Onarder  (Figs.  4a-4d)  and  -l-diffrac- 
tion-order  (Figs.  9a-9d)  meaaurementa  shows  that  much 
of  the  incident  energy  ia  re-emitted  into  the  diffraction 
order.  This  enhancement  in  the  dinraction-iader  inten¬ 
sity  is  mediated  hy  SPW  interactions. 

Finalfy,  the  diffraction-order  measurements  are  shown 
for  gratings  with  residual  SPW  coupling  effects  (Fig.  10a) 
and  without  SPW  effects  (Fig.  10b),  corresponding  to  the 
0-order  scans  of  Fig.  6.  For  the  shallower  grating,  there 
ia  a  decrease  in  intensity,  and  two  broad  minima,  at  18.2* 
and  24*,  are  observed,  with  the  minimum  at  18*  corre¬ 
sponding  to  the  maximum  in  the  0  order  (Fig.  6a)  and  the 
minimum  at  24*  corresponding  to  neither  the  minimum 
nor  the  maximum  in  the  0  order  (Fig.  10a).  For  the  deep¬ 
est  grating  measured,  there  are  no  residua)  SPW  effects 
but  rather  a  substantial  decrease  in  the  diffraction-order 
intensity  for  all  angles. 


1  1  1  »_i  4  2  0  K  I  '  o  .  5  i_i 


Fig.  7.  Cross-MCtion  SEM  pbotognphs  of  the  gratingi  used  for 
the  measuremenU  in  Fig.  6. 

either  the  nriinima  or  the  maxima  of  the  line  8hi4)e8  first 
dckreases  (Figs.  8a-8c)  and  then  increases  (Fig.  8d). 

Fot  deeper  gratings,  corresponding  to  the  0-order  scans 
of  Figs.  4a-4d,  the  diffraction-order  scans  are  shown 
in  Figs.  9a-9d,  where,  as  ak.ve.  Figs.  9e-9h  are  the  result 
of  theoretical  modeling  that  will  be  discussed  in  Sec¬ 
tion  8.  The  following  featmvs  can  be  noted  from  these 
measurements: 

(1)  There  is  an  increase  in  the  angular  separation  of  the 
two  minima  with  increasing  grating  depth,  in  contrast 
with  the  behavior  in  the  0-order  measurements,  where  the 
minimit  mow  toward  each  other  and  merge, 

(2)  There  is  a  decrease  in  the  angular  separation  of  the 
two  maxima  (Figs.  9e-9c)  and  an  emergence  of  a  single 
maximum  (Fig.  9d), 


5.  ENERGY-SUM  MEASUREMENTS 

For  these  measurements,  only  the  0  and  -1  orders  are 
propagating  within  the  angular  range  of  SPW  inter¬ 
actions.  Thus  we  can  evaluate  the  total  energy  coupled  to 


INCIDENT  ANGLE  (deg) 

Fig.  8.  -l-difTraction-order  scans  for  the  gratings  used  in 
Fig.  3.  The  left-hand  ootumn  presents  experimental  results;  the 
ri|^)t-hand  odumn  presents  theoretical  modding.  See  the  text 
for  details. 
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Fig.  a  -l-difTraction-order  acana  for  the  gratinga  in  Fig.  4. 
Tte  left'hand  column  preaenta  experimental  reaulta;  the  right* 
band  column  preaenta  theoretical  modeling.  See  the  text  for 
detaila. 


SPW’s  by  adding  the  energies  in  these  orders.  Fig* 
urea  lla-lld  show  these  energy  sums  for  the  measure* 
ments  of  Fig.  3  (0  order)  and  Fig.  8  (-1  diffraction  order). 
The  measurements  for  the  deeper  gratings  (Figs.  4  and  9) 
are  given  in  Figs,  lle-llh.  Note  that 

(1)  The  coupling  efficiencies  are  reduced  in  comparison 
with  the  0  order,  since  siune  of  the  incident  energy  is  also 
coupled  into  the  diffraction  order, 

(2)  The  coupling  of  energy  to  SPW’s  in  the  -2  order  is 
manifested  by  the  resonance  line  shape,  in  comparison 
with  the  les8-well*defined  line  shape  of  Figs.  3a-3^ 

(3)  There  is  a  gradual  increase  in  cou|ding  efficiency 
and  a  broadening  of  the  flrstKsrder  resonance  line  sh^fie 
followed  by  a  decreasing  efficiency  at  still  deeper  gratings, 

(4)  ‘ntere  is  a  decreasing  angular  separation  (Figs.  11a- 
llf)  followed  by  an  increasing  separation  for  deeper 
gratings. 

Both  the  -l*order  measurements  and  the  energy*sum 
measurements  show  a  well-defined  momentum  gq>;  Le., 
they  show  two  resonances  as  a  functkm  of  6  at  fixed  m 
that  first  i^iproach  each  other  ai^  then  diverge  as  the 
grating  depth  is  increased. 

The  experiments  exhibited  strong  dependence  on  the 
details  ci  the  grating  profiles.  For  some  cases,  with  ap¬ 
proximately  sinusoidal  jeoflles,  saddle-point  line  shapes 


were  not  observed  in  the  0-order  reflectances.  Also,  the 
diffraction-order  line  shapes  were  similar  to  the  0-order 
line  shapes  at  SPW  resonance  angles.  Not  unexpectedly, 
the  higher-order  Fourier  components  of  the  grating  pro- 
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Fig.  10.  —  l-diffraction-order  acana  for  the  gratinaa  uaed  in 
Pig.  6. 
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Fig.  11.  Energy-sum  (0-order  reflectance  plus  -l-diffraction- 
order)  scans  for  the  gratings  used  in  Figs.  3, 4,  8,  and  ft 
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Fig.  12.  O-order  renecUnco  and  -l-diffractioii'order  acana 
for  the  IKl-iini'deap  gratinga  (^iga.  4a,  9a,  and  lie)  as  tha  wav«> 
length  is  increased  Cram  686  to  607  nm.  The  teCt-haod  cidunin 
presents  the  0-order  measurements;  the  ri^ht-haad  oohunn  pre¬ 
sents  the  -l-dif0«ction-order  measurements. 

files  substantially  affect  the  first-  and  second-order  cou¬ 
pling  strengths  and  line  shapes. 

&  VARIABLE-mVELENGTH 
MEASUREMENTS 

For  the  95-nm-deep  grating  (Figs.  4a,  9a,  and  lie),  a  series 
c(  0-order  reflectance  and  -l-diffraction-order  measure¬ 
ments  was  carried  out  as  a  function  of  wavelength.  The 
output  from  a  dye  laser  was  tuned  through  the  580- 
633-nm  range,  sufficient  to  probe  the  (4-1,  -2)  minigap 
regioii. 

Figure  12  abows  the  0-order  reflectance  (Figs.  12a-12d) 
and  -l-diffractkm-order  (Figs.  12e-12h)  scans  as  wave¬ 
length  is  increased  from  686  to  607  nm;  the  vertical  scale 
for  the  diffraction-adtf  measurements  is  ten  times  more 
sensitive  than  the  corresponding  scale  for  the  0-order 
measurements.  Note  the  fdkming  essential  features: 

(1)  The  angular  separation  decreases  between  the  n  » 
■f  1  and  n  »  -2  coii^ingB,  with  the  angular  shifts  in  the 
second  ordw  being  almost  twice  those  in  the  first  m-der, 

(2)  A  wril-deflned  resonance  is  absent  in  the  0-<»rder 
seansoorreqMmdingtothea  <■  -2  coupling,  and  the  broad 
SPW  resonance  corresponding  to  the  n  «  4-1  coupling 
moves  toward  smaller  angular  values  (Figs.  12a-12<^ 

(3)  The  intensity  in  the  -1  diffraction  <wder  is  gnahially 
reduced  in  the  angular  region  between  the  two  couplings. 


(4)  As  in  the  previous  measurements,  the  -l-diffraction- 
ordCT  measurements  show  eomparaUe,  mirror-image  cou- 
jdings  for  the  two  SPW  resonances,  while  the  -2-order 
cou{ding  is  barely  observaUe  in  the  0-order  reflection. 

A  further  increase  in  wavelength  brings  the  two  cou¬ 
plings  stiH  closer,  and  scnne  interference  effects  are  ob¬ 
served  as  shown  in  Fig.  13.  The  major  features  are  the 
fcdknving: 

(1)  The  monotonic  shift  to  smaller  angular  values  in  the 
0-order  reflectance  scans  corresponding  to  the  n  »  4-1 
coupling  as  the  wavelength  is  changed  (Figs.  13a-13d), 

(2)  The  appearance  of  a  sharp  feature  at  23.8*  in  the 
-l-diffraction-order  measurements  (Figs.  13e-13f), 
which  a^reapooda  approximately  to  the  horizon  for  the 
appearance  of  the  -2  diffraction  order, 

(3)  The  appearance  of  a  slight  discontinuity  in  the  -1- 
diffractkm-order  line  shspe  (Fig.  13g), 

(4)  The  decrease  in  the  an^lar  separation  between  the 
two  coigdings  (Figs.  13e-13g), 

(5)  The  sharp  increase  in  the  diffraction-order  intensity 
as  the  two  couplings  cross  each  other. 

In  Fig.  14  these  measurements  are  extended  to  still 
longer  wavelengths,  beyond  the  (4-1,  —2)  minig^  r^on. 
The  m^jor  features  include 

(1)  The  appearance  of  a  =  -2  coiq>ling  in  the  O-order 
scans  (Figs.  14a-14d), 
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Fig.  IS.  Same  as  Fig.  12  hut  for  an  inwease  in  wavekngth  from 
610  to  617  nm. 
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Fig.  14.  Same  as  Fig.  12  but  for  an  increase  in  wavelength  from 
619  to  633  nm. 


(2)  The  s^ipearance  of  a  sharp  feature  at  24.2’  (Fig.  14e) 
corresponding  to  similar  features  in  Figs.  13e-13f, 

(3)  The  increase  in  angular  separation  between  the  n  = 
4-1  and  n  s  -2  couplings, 

(4)  The  increase  in  -1-diffraction-order  intensity  in  the 
region  between  SPW  couplings, 

(6)  The  reversal  of  the  resonance  line  shapes  in  the 
din'raction-order  measurements,  evident  from  a  compari¬ 
son  of  Figs.  12e-12g  and  14e-14g. 

Thus,  as  couplings  to  SPW’s  in  the  first  and  second 
orders  become  degenerate,  interference  effects  result  in 
an  enhancement  in  the  diffraction-order  intensity  over  a 
narrow  angular  range.  When  the  modes  cross  ea(±  other, 
the  resonances  observed  in  0  mder  show  a  sin^  minimum 
and  do  not  exhibit  a  momentum  gq>  (Fig.  16a),  whOe  a 
clearly  defined  forbidden  momentum  gap  is  observed  fw 
the  -l-diffraction-order  scans  (Fig.  16b).  The  magnitude 
of  this  gap  is  ~0.3*  and  is  in  good  qualitative  agreement 
with  the  momentum  gape  observed  ^  Heitman  et  al.** 

7.  THEORY 

A  variety  of  mathematical  formalisms  have  been  apfdied  to 
describe  cqrtical  interactions  with  poiodk  surfaces.  Inte¬ 
gral  meth^  were  first  proposed  by  Petit  and  Cadilhac,* 
Wirgin,”  and  Uretski.**  A  number  of  integral  fmmuhi- 
tkms  hiwe  been  investigated,  including  a  Oreen’s-functkm 
formalism  based  on  an  extinction  thecHrem  mechanism 
developed  fay  Glass  et  of.,**  Lake  et  and  Utagawa." 


In  addition,  differential  formalisms  devel(q;>ed  by  Nevi^ 
et  al.**  and  Numata**  have  also  been  used  to  study  SFW 
dispersion  relations.  All  these  methods  are  rigorous,  hcdd 
f(Hr  arbitrary  grating  profiles,  and  require  extensive  com- 
putation. 

Simpler,  approximate  formulations  have  been  intro¬ 
duced  starting  from  the  Rayleigh  hypothesis^  and  from 
the  extinction  theorem  formalism.**  These  were  demon¬ 
strated  to  offer  a  good  approximation  of  the  exact  calcula¬ 
tion**  and  to  Hrst-order  coupling  experiments.*  While 
these  formulations  are  lacking  in  ultimate  rigor  and  accu¬ 
racy,  they  ofler  a  good  intuitive  physical  picture  of  the 
interaction  with  a  minimum  of  computational  difficulty. 
Here,  the  extension  of  this  simple  formulation  for  the 
minigap  region  is  briefly  described.  More  details  are  pre¬ 
sented  in  Appendix  A. 

For  TM-polarized  light  incident  on  a  sinusoidal  surface 
with  an  excursion  of  2u,  a  strai^tforward  but  involved 
algebraic  procedure  results  in  the  coupled  linear  differen¬ 
tial  equations 


Jp(k.u)  +  2Bl{-iy-Jp-.(anU)  -  2BTJp.,(B,u)  =  0, 

A  m 

(3) 

(a,  -  pg^y,{k,u) 

-  2  [a.  -  (P  -  n)g^jB;(-l)'-V,-.(o.u) 

-  ~  2  P.-ip-  rr)g^jB:J,-.O.u)  =  0,  (4)  0 
where  all  variables  are  defined  in  Appendix  A. 


SPW  COUPUNG  IN  THE 
(-H.-2)  MINI-GAP  REGION 


INCIDENT  ANGLE  (deg) 


Fig.  1&.  Plot  of  the  resonaiice  bdiavior  observed  in  Figs.  12-14. 
The  O^irder  reflectivity  (a)  abowe  a  tingle  dominant  resonance 
becoming  a  single  bro^  minimum  in  ^  minigap  region.  In 
ewtraat,  the  -l-diffractkm-crder  scans  (b)  show  a  weU-deflned 
momentum  gap  corresponding  to  the  (41,  -2)  minigap  region. 
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Fi*.  16.  Comparison  between  theory  and  experiment  for  a 
grating  with  a  nearly  sinusoidal  profile.  The  left-hand  column 
presents  experiment^  results;  the  right-hand  column  presents 
theoretical  modeling. 

In  the  derivation  of  Eqs.  (3)  and  (4)  no  small-signal 
approximation  is  made  for  either  the  grating  depth  m 
t^  diffraction-order  amfditudes.  From  the  dispersion- 
relation  curve  in  Fig.  1,  for  the  present  case  the  o^  reso¬ 
nant  terms  in  the  above  expansions  [Eqs.  (3)  and  (4)]  are 
n  =  +1  and  n  =  -2.  The  »  =  -1  term  is  Ae  permitted 
diffraction  order  and  is  strongly  influenced  hy  SFW  inter¬ 
actions.  The  uiHnite  series  was  truncated  at  n  =  0,  -1, 
1,  -2,  which  resulted  in  an  8  x  8  matrix  and  considerable 
simplincation.  Previous  analyses  included  many  more 
terms  and  provided  signiHcantty  greater  accuracy.  Since 
the  experimental  profiles  were  not  precisely  sinusoidal, 
additional  Fourier  conqwnents  were  included  in  Ae  grat¬ 
ing  profile.  It  was  shown  by  Van  den  Berg  and  Fokkema* 
that  the  Rayleigh  IiypoAesis  would  still  be  valid  if  a  pro¬ 
file  could  be  describe  by  a  finite  Fourier  series.  The 
mathematical  analysis,  however,  becomes  progressively 
more  complex  as  the  number  of  higher-order  Fourier  com- 
pmients  is  increased.  For  the  present  purpose  the  surface 
profile  is  redefined  to  include  only  on«  ad^tumal  Fourier 
component: 

X  ■  Hi  sinfgy)  -I-  u»  ria(2gy).  (5) 

In  general,  uj  is  much  smaller  than  ui.  For  this  {xofile 
we  obtain 

+  i)|A(As)  “  0  (*) 


D(P)  +  2  ^(P>  n)  -  —  2  f(P. «) 

•  ^  m 

+  sfolp.*)  +  SAp.".*)  -  —  =  0,  (7) 

.-iL  •  «  •  J 

where  A,  B,  C,  D,  E,  F,  G,  J,  and  L  are  defined  in 
Ai^ndix  A. 

Equations  (6)  and  (7)  reduce  to  Ecp.  (3)  and  (4)  for  u*  = 
0.  The  terms  in  these  equations  are  of  the  general  form 
J.J,  -,  in  the  numerical  calculations  terms  are  retained  iq> 
to  s  s'  =  3;  i.e.,  c/o</o,  <focls,  end 

B.  COMPARISON  BETWEEN  THEORY 
AND  EXPERIMENT 

In  a  previous  analysis  of  first-order  coupling  to  SPWs  it 
was  shown  that  Ae  best  Ht  to  Ae  experimental  results  was 
obtained  for  a  Ag-fllm  dielectric  constant  of  (-13.65;  0.82) 
(Ref.  5);  this  value  was  used  without  furtl»r  a4iustment. 
In  Fig.  16  experimental  and  theoretical  reflectance  scans 
versus  angle  are  shown  for  an  854-nm-period  grating  wi  A 
a  nearly  sinusoidal  profile.  The  Aeoretical  calculations 
were  carried  out  wiA  the  use  of  Eqs.  (6)  and  (7).  Note 
that  Ae  calculated  -f  1-order  coupling  resonance  line  shape 
in  Ae  -l-diffraction-mrder  scan  (Fig.  16b)  is  similar  to 
the  0-order  line  shiq>e  (Fig.  16a)  and  that  the  two  reso¬ 
nances  are  conq>arable,  although  only  >  smaD  enhance¬ 
ment  rather  than  a  coupling  dip  at  Ae  -2-order  resonance 
is  observed  in  Ae  0-order  scan.  The  Aeoretical  line 
shapes  (Figs.  16c  and  16d)  closely  follow  the  experimental 
line  shapes,  although  Ae  observed  resonance  widths  are 
slightly  larger  than  those  predicted. 

For  grating  profiles  deviating  from  a  sinuaoidal  shape, 
U|  and  Ut  were  adjusted  to  produce  line  shapes  closest  to 
those  experimentsdly  observed.  In  Fig.  17  the  experimen- 
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Fi».  17.  C(nnpariaoa  between  Aeoty  and  experinent  (tar  a  pro¬ 
file  showing  a  ali|d>t  deviation  fkmn  a  sinusoidal  shape.  The  left- 
hand  column  pwaents  experimental  results;  the  right-hand 
column  presents  theoretical  modriing. 
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tal  and  calculated  plots  for  a  grating  with  a  765-nm  period 
are  shown.  This  grating  deviated  from  a  sinusoidsJ  pro¬ 
file  and  showed  appreciable  coupling  in  the  second  order 
(Fig.  17a).  In  the  -1  diffraction  order  sharp  spikes  in  in¬ 
tensity  at  SI’W  resonance  minima  are  observed  (Fig.  17b). 
The  calculated  0-order  plot  (Fig.  17c)  is  similar  to  the  ex¬ 
perimental  plot  (Fig.  17a).  The  calculated  -1-diffraction- 
order  plot  (Fig.  I7d),  although  showing  intensity  spikes  at 
both  SPW  resonances,  fails  to  display  the  experimentally 
observed  dispersionlike  minima,  shows  unequal  intensi¬ 
ties  at  the  SPW  resonances,  and  does  not  replicate  the  ob¬ 
served  diffraction  efficiency  enhancement  in  the  region 
between  the  SPW  resonances. 

Thus  this  simple  model  includes  many  of  the  m^jor 
features  observed,  but  for  a  complete  description  higher 
Fourier  coefficients  must  be  included,  as  we  can  see 
in  comparing  the  experimental  and  theoretical  plots  in 
Figs.  3a-3d  and  3e-3h  above.  Note  the  following  features: 

(1)  The  resonance  line  shapes  in  the  second  order  are 
different  (Figs.  3a-3c  and  3e-3g), 

(2)  The  observed  grating  dept^  are  much  larger  than 
those  predicted  (Figs.  3a-3c  and  3e-3g), 

(3)  The  calculated  angular  shifts  are  much  smaller  than 
those  observed. 

The  large  difference  in  grating  depths  and  the  second- 
order  line  shapes  is  due  to  the  omission  of  higher  Fourier 
coefTicients  in  the  grating  profile  that  result  in  cancella¬ 
tion  effects  in  second-order  coupling.  Specifically,  the 
model  at  present  does  not  include  a  third-order  coeffi¬ 
cient.  Yet,  as  we  discussed  above,  it  is  just  this  Fourier 
component  of  the  grating  that  resonantly  couples  the 
counterpropagating  SPW  modes.  For  deeper  gratings 
the  comparison  between  observed  (Figs.  4a-4d)  and  pre¬ 
dicted  (Figs.  4e-4h)  reflectance  scans  show  the  foUowing 
features: 

(1)  There  is  an  evolution  to  a  broad  minimum  (Figs.  4e- 
4h), 

(2)  The  calculated  and  observed  depths  are  similar 
(Figs.  4a-4c  and  4e-4g), 

(3)  The  calculated  angular  shifts  are  smaller  than  those 
observed, 

(4)  The  calculated  energy  coupled  to  the  diffraction 
order  is  much  higher  than  that  observed. 

Note  that  the  calculated  plot  in  Fig.  4h  is  fw  a  sinusoidal 
profile,  while  the  actual  profile  shows  considendde  devia¬ 
tion  frinn  a  sinusoidal  shape  (Fig.  6d).  These  conqMri- 
acms  emphasize  the  impwtance  of  knowledge  of  the  exact 
profile  fw  calculation  of  its  reflectance  behavior. 

The  comparison  of  calculated  (Figs.  8a-8d)  and  ob¬ 
served  (Figs.  8e-%)  diffraction-order  line  shapes  shows 
the  folkranng  features: 

(1)  There  is  a  lack  of  energy  enhancement  in  the  region 
between  intensity  spikes  (Figs.  8e-8h), 

(2)  The  calculated  spikes  in  intensity  are  sUgfatly  asym¬ 
metric  (Figs.  8e-8g), 

(3)  The  calculated  energy  in  the  diffraction  ordo*  is 
hi^to’  than  that  observed. 

For  deepo-  gratings  a  similar  conq>ari8on  of  observed 
(Figs.  9a-9d)  and  calculated  (Figs.  9e-9h)  diffraction- 
ordu  line  shi4)es  shows  the  following  features; 


(1)  An  energy  enhancement  for  the  calculated  line  shapes 
in  the  region  between  SPW  resonances  (Figs.  9e-9h), 

(2)  The  appearance  of  shallow  minima  (Figs.  9e-9h), 

(3)  The  appearance  of  an  asymmetric  line  shape  in 
Fig.  9g, 

(4)  The  sharp  decrease  in  signal  from  ~10*  to  15*  in 
Fig.«i. 

Thus  the  con4>ari8on  of  theory  and  experimental  results 
indicates  that,  although  the  theory  manages  to  show  rnsjor 
featiues  of  our  data,  many  Hner  details  are  not  observed. 
A  more  detailed  theoretic  investigation,  including  both 
hitler  Fourier  coeRlcients  and  a  more  rigorous  formal¬ 
ism,  is  necessary  to  provide  a  detailed  fit  between  theory 
and  experiment. 

9.  SUMMARY 

A  q^matic  experimental  stu(fy  of  interactions  between 
SPW  couplings  in  the  (+1,-2)  minigap  r^ons  is  reported. 
The  experiment  covers  the  entire  range  of  SPW  inter¬ 
actions,  ranging  from  weak-coupling  to  stronger-coupling 
effects  and  to  deeper  profiles  showing  absorption  rather 
than  SPW  effects.  The  simple  theoretical  model,  al- 
tbou^  miasing  many  of  the  finer  details,  shows  the  nugor 
features  of  oiir  experimental  data.  One  can  develop  a 
better  model  by  including  several  Fourier  coefTicients  in 
the  grating  profile. 

Several  extensions  of  this  work  are  immediatebr  appar¬ 
ent.  For  larger  grating  periods  there  are  more  diffrac¬ 
tion  wders,  and  it  would  interesting  to  see  how  energy 
is  distributed  in  these  orders  at  SPW  resonance  angles. 
Also,  the  appearance  of  momentum  gaps  in  the  0  order 
depends  on  the  strength  of  interactions. 

As  we  noted  above,  there  has  been  much  discussion  of 
tbe  physical  significance  of  momentum  aixl  energy  gaps  in 
the  minigap  regions.  As  this  stuc^y,  as  well  as  others,*"^* 
demonstrated,  the  respcmse  surface  is  quite  conplex  in  the 
vicinity  of  a  minigap,  and  a  complete  analysis  of  the  ex¬ 
perimental  situation  is  necessary.  Most  theoretical  treat¬ 
ments  of  SPW  effects  on  gratings  can  be  cast  as  a  linear 
algebraic  proUem  of  tbe  general  form  ME  =  8,  where  the 
vector  E  represents  the  response  fields,  M  is  the  medium 
response  matrix,  and  8  is  the  source  of  excitation  vector. 
Formally,  this  is  solved  as  E  =  M~'8,  where  each  term  in 
the  inverse  matrix  is  proportional  to  (det  M)~‘.  In  many 
situations  knowledge  of  this  determinant,  and  hence 
of  the  diqierskm  rdatkm,  is  sufficient  for  an  undostand- 
ing  of  the  main  features  of  the  respmise.  For  example,  for 
the  case  of  a  sin^e  SPW  resonance,  the  experimental  dis¬ 
persion  relation  obtained  from  the  resonance  angles  and 
the  linewidths  of  the  0-order  coupling  dips  are  in  reason¬ 
able  agreement  with  that  calculate  from  the  determinant 
of  M;  higher-order  terms  can  be  added  to  evaluate  the  ef¬ 
fect  of  deeper  gratings.*  In  the  minig^>  regi<m  this  pro¬ 
cedure  clearly  is  not  apfsof^te;  for  exanple,  neither  of 
the  extrema  observed  in  tbe  first-order  coufding  coincides 
with  tbe  0-order  coupling  resonance  angles.  In  this  situ¬ 
ation  it  is  not  possible  simply  to  extract  a  dispersion 
relation  from  the  experimental^  observed  line  shapes.  A 
full  treatment  of  the  qrstem  response  is  necessary  fw  in¬ 
terpretation  of  the  measurements  in  terms  of  modes  and 
diqiersion  relations. 
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APPENDIX  A 

For  a  TM-polarized  plane  wave  incident  (m  a  einuaoidal 
interface  we  can  deHne  a  (magnetic)  B  field  below  and 
above  the  grating  I9  using  tlw  Rayleigh  hypothesis: 

»  exp(»(*,y  +  *.*)]  +  ~  «**)! 

m 

Z  <  0, 

B"(>iz)  =  S  Sr  exp[i(ir,y  +  /3,z)],  z  >  0, 

m 

where  A,  =  Ao  sin  8,  k,  =>  kt  cos  Aq  »  ai/c,  k^  ^  ky  +  ng 
(n  «  0,  ±1,  ±2),  g  •=  2v/d  (d  is  the  grating  period),  and 

a,  =  i(A.*  -  Ao*)*^.  B.  “  «•(*.*  -  cV)*^. 

Fw  a  sinusoidal  surface  the  profile  is  given  by 

f{y)  -  u  ainigy), 

so  that,  by  applying  the  generating  function  for  Bessel 
functions,  we  obtain 

«tp(»yAy)l  ®  exp(iy«  sin(gy)]  =  2  “I>(wy)«/,(ya). 

By  applying  the  boundary  conditions 
B*(y,z)  =  B-(y,z)l.-A,„ 


P(P.»)  -  B^B.  -  (P  -  /i)g(A./B.)l 

<f»-a(B«Kl)«fo(B«U|)  -  gA.UiBr 
[•(»-i«-»(B»Bi)  +  <f^-a+l(BiiUl)]c/»(B«>tt)i 
0(A*)  -  t*.  -  (P  -  2*)g(V*.)K-i.(*.«i)«/.(*.«,) 

+  (-m*.  -  (p  +  2»)g(ky/k.)\J,,u(k.u,) 

X  t/,(A,Kt)  ~  8kjUiJf,u-ii.k,Ui) 

+  ('llV/fafa-iiAiUilJc/tiA.ua)  -  gAyUi 
^  {.’fp-tM*iik,Ui)  +  (— l)*€^.tjj-fj(Aji<j)]«/,(A,Kj), 

/(p,n,s)  =  -(-1)'”*B5«.  -  (p  -  n  -  28)g(A,/o,)] 

X  J,-.-n{o.Ui)«f.(a.«i)  +  (-1)*BI 
^  [a.  -  (p  -  n  +  2s)g(A,/o,)] 

^  •f»-»+**(o«ai)tf,(a,Ua)  —  gk„UiB^ 

X  [«/p-»-*,-*(«,Ui)  +  (-VfJf 

-w+a*-  j(o.ltl)] 

X  <7.(a.Ut)  -  gA,«,B:[J,-,.j,*»{o«Ki) 

+  ('t'l)^e^-a42t'<'2(0||U|))c/«(a|U}) , 

Hp,n,»)  =  BrtBa  -  (p  -  n  -  2s)g(A,/B,)] 

X  Jp-.-iAfi.Ur)J.(fi,Ut)  +  (-D-B: 

[B»  -  (P  -  n  +  2s)g(A./B.)I 
X  ‘fp-»+».(B»“i)*^t(B.“*)  -  g*«u»Br 
^  [«fp-«-t*-i(B»“i)  +  -»4*f-2(B*ai)l 

^  •f«(B»tti)  ~  g^»aiB]J(e7p-i,_*j4i(B»ai) 

+  (~l)^*^i-»-H*-n(B»“l)]«f»(B»M»)- 


where  9/dn  is  the  normal  derivative  to  the  surface  and  is 
defined  by 


dn 


we  arrive  at  Eqs.  (3)  and  (4)  above.  For  a  more  complex 
sinusoidal  surface  the  profile  is  given  by 


Ay) 

and  we  obtain 


«i  ein(gy)  +  ut  sin(2gy). 


MrfiyAy)]  “  2  «pI»(p  +  29)ByK(y«i)«f»(yMi). 
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Again,  applying  the  boundary  conditions  leads  to  Eqp.  (6) 
and  (7)  abo^  where  the  constants  are  given  by 

Mp,»)  -  [Jp-u(k.ux)  +  {-irJp*u{k.Ux)}J.{k^,). 
B{p,n,t)  =■  (-ir*BJJ,...fc(a.u,) 

+  (“i)^*f»-»+i»(«ii«i)W(««*<i)i 

eXAn.s)  »  Br[c;,...«.(B.Ui) 

+  (-irJ,-.4*.(B.Mi)W(B.B*). 

Dip)  -  [Jf,  -  pg(A,/A.)]J,(A.Ui)Xo(A.»i)  -  gkjUt 

Eip,n)  -  -i-lf-Bia,  -  (p  -  n)g(A^0] 

X  J^.,ia,Ut)JUa,u^  -  gk,u,Bl 
X  [•7,...i<a.Bi)  +  Jp.,*dei,Ui)}Jda,Ut), 
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